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Abstract. The concentration and distribution of quantum entanglement is an 
essential ingredient in emerging quantum information technologies. Much theoretical 
and experimental effort has been expended in understanding how to distribute 
' entanglement in one-dimensional networks. However, as experimental techniques in 

quantum communication develop, protocols for multi-dimensional systems become 
essential. Here, we focus on recent theoretical developments in protocols for 
distributing entanglement in regular and complex networks, with particular attention 
to percolation theory and network-based error correction. 
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1. Introduction 

The idea of quantum entanglement has a long history [H El El S], although an 
intensive search for a comprehensive theory of entanglement only arose with quantum 
information theory. This search grew out of the realization that quantum entanglement 
is an essential resource for developing information technologies that are radically 
different than those possible in a purely classical world. In fact, when two physical 
systems are sufficiently entangled, they exhibit correlations that are stronger than 
possible with any classical theory. These strong correlations can then be exploited by 
cryptographic, communication, or computation protocols. As with classical information 
theory, there is a fundamental need to understand how to distribute information, that 
is, how to transmit a signal between two parties. But quantum mechanics imposes 
severe limitations, both fundamental and practical, on copying, encoding, and reading 
information. Thus, distributing quantum entanglement is an extremely challenging 
problem. This review presents work that attempts to meet that challenge. 

Knowledge of the network in which this entanglement distribution will take place 
is still in a nascent stage. Naturally, most attention was initially focused on one- 
dimensional setups [5l El [7]. However, it is natural to consider distribution on multi- 
dimensional or otherwise more highly-connected networks, whose structure may be 
designed explicitly for distribution or may imposed by geographical constraints. Obvious 
examples of connections to existing bodies of work that will become increasingly 
important are the classical theory of complex networks, particularly those with an 
internet-like structure jH [91 [101 [11] . From another point of view, the production and 
manipulation of entanglement on micro- or nano- scale networks using, for instance, 
photonic crystals is progressing rapidly [121 H3] • In this case, we will likely be presented 
with regular arrays of elements that can be entangled. Both of these kinds of systems, 
and others not yet imagined, will require an understanding of entanglement distribution. 

In any communication application, information is encoded in the state of a physical 
system. But encoding information in quantum mechanical states poses new challenges. 
For instance, communication elements, such as repeaters, amplifiers and memories, rely 
on repeatedly measuring a system. The difficulty becomes immediately apparent when 
we recall one of the most fundamental features of quantum mechanics: in general, a 
measurement reveals only partial information on the state of the system, while at the 
same time suddenly "collapsing" it into a new state in which some or all of previous 
information is lost. On the other hand, quantum mechanical entanglement opens new 
possibilities in manipulation of information that are fundamentally impossible in the 
classical theory. 

Before giving a detailed description of entanglement, we briefly review a few of the 
most important applications of this phenomenon; we refer the reader to [H] for a more 
complete review. These examples will repeatedly refer to the fundamental quantum 
system for quantum information science, the qubit, which is used as a quantum analogue 
of the classical binary digit or bit. The qubit is an abstraction that is realized by any 
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two-level quantum system: the spin of an electron or neutron, the polarization of a 
photon, or the energy of an atomic electron in a resonant field, just to name a few. 
These systems have a single, two-dimensional degree of freedom, which means that, for 
a given orientation of the measurement device, a measurement always gives one of two 
results. 

Quantum teleportation It has been proven that the unitary evolution of quantum 
mechanics implies that an unknown quantum state cannot be duplicated or cloned [15J. 
However, an unknown quantum state can be transported over an arbitrarily long 
distance as long as an auxiliary perfectly entangled pair of particles (also called a Bell 
pair) and a classical communication channel are established over the same distance [16J. 
The entangled pair is created via a local temporary interaction between two qubits, 
which are then separated by the desired distance. The procedure is as follows. Two 
distant parties, traditionally called Alice and Bob, share an entangled pair of qubits, 
while Alice has an additional "data" qubit that she wishes to send to Bob. Alice 
performs a certain measurement on the two qubits she possesses and communicates 
the result classically to Bob. Bob then applies a transformation on his qubit in a 
manner prescribed by the message from Alice. The result is that Bob's qubit is now in 
the original state of the data qubit of Alice, which meanwhile has lost its information 
content. 

Quantum distributed computing Distributed computing consists of several nodes that 
do computations independently while periodically sharing results. Entanglement 
appears in several places in quantum distributed computing protocols, including in the 
input states or in the communication channels used in sharing results between nodes. 
Attempts have been made to design quantum distributed computers so that limited 
entanglement resources are spread in an optimal way among components [T7] . 

Quantum key distribution Classical public key cryptography schemes are widely used, 
for instance, in internet security algorithms. Entangled pairs can be used to securely 
generate and distribute the classical private key necessary for these schemes [T51 EES]- 
Using a quantum protocol, Alice and Bob generate a series of random bits, sharing 
knowledge of the results, but preventing others from eavesdropping. To this end, 
they initially share a number of Bell pairs, which they measure sequentially, each time 
choosing an orientation randomly from a pre-determined set. They can use a portion 
of their results to compute a bound on the amount of eavesdropping (or noise) that has 
taken place, and another portion to generate the key. 

Superdense coding If two parties share a Bell pair, then two classical bits of information 
can be sent from one party to the other one, even though each party physically possesses 
only one qubit [20]. To accomplish this, Alice applies one of four previously agreed- 
upon unitary operations to her qubit. An unitary operation, in contrast to measurement, 
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transforms the qubit in a non-destructive and coherent way. This transforms the Bell 
state to one of four orthogonal states that together form the Bell basis. Alice sends her 
qubit to Bob, who then performs a joint measurement on both qubits, thus distinguishing 
reliably among the four messages Alice can send. 

Each of these applications requires entangled pairs of particles to be generated and 
distributed between two distant parties. Currently, the main technological difficulty is to 
create remote entanglement which, in most experiments, is achieved by sending polarized 
photons through optical fibres. In fact, due to noise, scattering, and absorption, 
the probability that the quantum information contained in such photons reaches its 
destination decreases exponentially with the distance. Another experimental challenge 
is to transfer the quantum state of the photon onto that of a quantum memory, such 
that the entanglement can be manipulated and stored; see [21] for a recent review on 
quantum memories. 

Outline We first review some well-known results on entanglement and describe the 
operations that allow one to propagate quantum information in a network, such as 
entanglement distillation and entanglement swapping; see section 12.1.21 and section 12.21 
Based on these operations, the quantum repeaters enable entanglement to be generated 
over a large distance in one- dimensional networks; see section 12.31 However, in order to 
obtain reasonable communication rates, they require an amount of physical resources 
that increases with the distance, which is out of reach with present technologies. Since 
quantum communication networks will typically be spread over the Earth's surface, one 
can ask whether the higher connectivity of the stations or nodes of the network will 
be of some advantage over a one-dimensional setup when distributing entanglement; 
see section [3J The first protocol exploiting this fact uses ideas of percolation in two- 
dimensional graphs: for pure states, if enough entanglement is generated between 
neighbouring stations, then it can be propagated over infinite distance; see section 13.21 
However, for this result to hold, quantum operations have to be applied perfectly on 
the system. Obviously, this cannot be achieved in practice; but for some specific types 
of noise, entanglement percolation still allows entanglement to be generated between 
infinitely distant stations; see section 13.31 

We will see that all percolation strategies need, in the end, perfect operations to 
be applied on the system. Their study is certainly useful; for instance, in establishing 
fundamental limits. Still, we must ask if the higher connectivity of realistic quantum 
networks will be of practical use. This is indeed the case, but radically different 
protocols, based on error correction, must be designed; see section HI In fact, while 
entanglement percolation relies on the existence of one path of perfectly entangled states 
between any two stations, network-based error correction extracts its information from 
all paths connecting the two stations. Contrary to the quantum repeaters, no quantum 
memory is needed. 

Initially applied to regular graphs, the study of entanglement distribution has since 
been extended to complex networks; see section This is a natural generalization 
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because present communication networks exhibit a complex structure. Furthermore, 
while mostly focused on pure states, recent work concerns the manipulation of 
entanglement in noisy quantum complex networks; see section 15.41 These results 
emphasize the fact that the quality of the quantum communication between two stations 
will depend greatly on our understanding of the interplay of the network topology and 
the quantum operations available on the system. 
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2. Concepts and methods 

In this section we describe some well-known results on entanglement and the basic 
operations that are used to propagate quantum information in a network; see [22] for 
a thorough treatment of this material in a pedagogical setting. A reader familiar with 
quantum information theory may skip this part and jump to the discussion of quantum 
repeaters in section 12.31 or percolation in section [3J 

2.1. Quantum states 

Here we give a very brief account of the concepts necessary to understand entanglement 
below. We denote the state of a single qubit by 



where |0) and |1) represent a choice of basis, called the computational basis, in a Hilbert 
spacjj]. We say that the system is in a coherent superposition of the two basis states. 
We can measure the state of the qubit in this basis, which corresponds, for instance, 
to the orientation of magnets acting on the spin of electrons in the laboratory. The 
probability that we measure is ao, while the probability that we measure 1 is a±. We 
must have a + «i = 1, as these are the only possible outcomes. Upon measurement 
in the computational basis, the state of the system collapses into one of the two basis 
states, say |0). If we now repeat the same measurement, we have a — 1 and a± = 0, 
so that we obtain |0) with probability 1. However, if we rotate our magnets to an 
orientation corresponding to a different basis and then measure, we must expand |0) 
in that new basis in order to calculate the probabilities of obtaining each of the two 
possible results. 

2.1.1. Measurements and quantum evolution While the study of quantum measure- 
ment is a broad and deep subject, here we only need to introduce a few ideas to discuss 
entanglement distribution. Measurements occur in the laboratory, but the computa- 
tional tool to predict their result associates with each type of measurement a set of 
linear operators on the state space. The measurements described above are called pro- 
jective measurements, which are defined by a collection of projectors Ej onto subspaces 
of the state space, each of which is associated with a possible measurement value. After 
measuring a value, we know that the system has collapsed into a state in the subspace 
corresponding to the associated projector. The projectors are orthogonal, and the re- 
quirement that we must get some result implies that their sum is the identity. The 
maximum number of orthogonal projectors is equal to the dimension of the state space 
and corresponds to a complete measurement. On the contrary, we use an incomplete 
measurement, with fewer projectors, when we only want some information from a state 
while preserving some property common to all subspaces. Finally, we will sometimes 

f We neglect here the relative phase in the superposition of (|2.ip which will be inessential for the full 
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make use of the generalized measurement in which the condition that the operators are 
orthogonal is relaxed. Measurements are then defined by a collection of semi-definite 
positive operators that sum to the identity [23] . 

The other component of quantum theory that transforms states is evolution. 
Evolution is also represented theoretically by operators, but in contrast with 
measurement, an initial state is transformed into a final state in a deterministic way. 
Operators representing evolution must have the algebraic property of unitarity in order 
to preserve probability density. These operators are commonly referred to as unitaries. 
Quantum information protocols are built from a combination of measurements and 
unitaries on quantum states. 

2.1.2. Entanglement To illustrate the basic idea of entanglement while making a 
connection with the remainder of this review, let us consider two qubits, i.e. two systems 
each of which consist of a two-level quantum state. Two qubits live in a four-dimensional 
state space whose computational basis is 1 00) , |01), |10), |11), where the first binary 
digit corresponds to qubit A and the second to qubit B. Mathematically, these basis 
vectors are tensor products of vectors in the local spaces: \ij) = \i) Cg> If the two 
systems have never interacted directly or indirectly, then each one has an independent 
description of its state, as in (12. ip . It is possible to find local bases for each qubit such 
that the joint state is one of the four computational basis states of the joint space. These 
states are called product states. Measurements and operations applied on one system 
have no effect on subsequent measurement outcomes on the other system. However, 
suppose an interaction between the systems is turned on, then off. In general, there is 
no local basis such that the state of the system after the interaction can be written as 
a factor of states of the subsystems. A pure bipartite quantum state is then said to be 
entangled if it cannot be written as a tensor product. That is, it cannot be written as 



which form a basis of two-qubit systems. These states are of central importance as they 
are maximally entangled states. The projective measurement corresponding to the four 
states of the Bell basis is called a Bell measurement. 

Just as two two-dimensional systems form a four-dimensional system, in general, 
the space describing a collection of quantum systems of dimensions n 1; n 2 , . . . , has 
dimension iii x 112 x . . . x %. However, generalizing the case of two qubits, we may 
partition any of these composite spaces into two subspaces, and viewed this way, the 
system is called a bi-partite system. The study of bi-partite entanglement, that is 
entanglement between the two subspaces, is far better understood than the more general 




(2.2) 




(2.3a) 



(2.36) 
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case of multi-partite entanglement, and we shall mostly be concerned in this review with 
bi-partite systems. 

Until now, we have discussed only what are known as pure quantum states, in which 
the probability of measuring a value is of local quantum mechanical origin. However, 
the generic case is that the system is entangled with other systems which we cannot 
measure. It turns out that to predict local measurements in this case, we can assume 
that the local state is effectively a classical distribution of several pure states. This 
classical ensemble of pure states is known as a mixed state. These states are described 
not by vectors in the Hilbert space of pure states, but rather by linear operators known as 
density matrices, which act on the same Hilbert space. For instance the density matrix 
of a pure state is given by the outer product of the pure state vector with itself, which is 
denoted by \<p)(4>\ ■ Rules of quantum mechanics together with classical statistics imply 
that density matrices must be positive and have trace one. 

The notion of entanglement is easily generalized to mixed states [21]. In this case, 
a bi-partite system is entangled if and only if its joint density matrix pab cannot be 
written 



that is, it is not a mixture of product states. 

Measures of entanglement If one qubit of a Bell pair is measured in the computational 
Z basis, that is, if it is projected onto the eigenstates |0) and |1) of the Pauli Z matrix, 
then the measurement outcome will be or 1 with equal probability. If the second qubit 
is then measured in the same basis, the outcome will be completely determined by the 
first result. But, if we instead were to measure each qubit in the same arbitrary rotated 
basis, the same correlation between outcomes would be seen. In fact, it can be shown 
that the Bell states are the maximally-correlated two-qubit states; they are maximally 
entangled. In operational terms this means that they can be used to teleport perfectly 
the state of exactly one qubit. In practice, however, the entanglement between two 
qubits is never perfect, so that partially entangled states have to be considered. In what 
follows we present some common measures of entanglement. 

In the pure-state formalism, local bases may be found so that any two-qubit state 
may be written 



where the two Schmidt coefficients (po and <pi satisfy ipo + (p\ = 1 and ipo > (p\ (by 
convention). If one of the coefficients (po or (pi vanishes, the system is separable. If 
V?o = fi — 1/2, it is maximally entangled. Otherwise, the system is said to be partially 
or weakly entangled. A measure of entanglement of this state is given by 




(2.4) 




(2.5) 



E(<p) = 2(p! e [0, 1] 



(2.6) 
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which corresponds to the optimum probability of successfully converting |<p) into a 
perfect Bell pair by local operations [22]; see section [2.2.21 Another common measure 
of entanglement, the concurrence [26], reads in the case of pure states: 

C(¥>) = 2VwIe[0,l]. (2.7) 

The effects of imperfections, or noise, on a realistic system are taken into account 
if one considers mixed states. In this case, the state of two qubits is given by a four-by- 
four density matrix, which requires, in general, fifteen parameters. Any density matrix 
of two qubits a can be transformed by local operations to the standard form 

p w (x)=x|$+>($ + | +^id 4 , (2.8) 

where id 4 is the corresponding identity matrix. This process is known as depolarization, 
and pw is called a Werner state. One sometimes explicitly expands this expression in 
the Bell basis B: 

. . f \-F \-F \-F\ . s 

MF)=(F,— — — ) b , (2.9) 

where the components correspond to the weight of each Bell state in the mixed state 
decomposition of pw(F). The fidelity of the state is not altered by the depolarization 
procedure: 

F(p w ) = ($ + | pw |$ + > = (3x + l)/4. (2.10) 

This fact is important because the concurrence of the Werner state pw is related to the 
fidelity via 

C(p w ) = max{0,2F(p w ) - 1}. (2.11) 

The state p w is entangled if and only if its concurrence is strictly positive, that is if 
x > 1/3; or equivalently, if 

F>F min = 1/2. (2.12) 

Note that ( 12. 8 p emphasizes the fact that a Werner state is a mixture of a perfect quantum 
connection and a completely unknown state. This observation is important, for instance, 
when discussing the limitations of entanglement propagation in noisy networks; see 
section 13.41 The Werner state is an example of a full-rank state, that is, one with no 
vanishing eigenvalues. Full-rank states are important because they model any state 
found in the laboratory, as they result from any general noise model. Moreover, as we 
shall see in the next section, these states are the most difficult to distill. For instance, 
a single copy of a pure (i.e., rank-one) state may be probabilistically converted to a 
Bell state. Likewise, a finite number of copies of a certain class of rank-two states may 
be also probabilistically converted. But an infinite number of full-rank states, Werner 
states for instance, are required to distill a Bell state. 

The discussion to this point may leave the reader with the impression that 
characterization of entanglement is a relatively uncomplicated task. In fact it is a 
difficult and deep question, especially for mixed states. On one hand a separable state 
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Figure 1: A quantum network. The stations (nodes) initially share some partially 
entangled pairs of qubits (links); for simplicity, we sometimes omit explicitly drawing 
the qubits at the stations. By LOCC, we mean that any quantum operation can be 
applied on the qubits of a station, but that only classical information is exchanged 
between the nodes. 



can be written as a classical mixture of projectors onto product states. On the other 
hand, there is no general algorithm to determine whether a given density matrix can be 
written in this form, although significant progress has been made [27| fT¥j . 

2.2. Entanglement manipulation in quantum networks 

In the introduction, we have seen that entanglement is a valuable resource for a variety 
of quantum information applications. It is thus essential to understand the non-trivial 
task of creating and distributing entanglement between distant parties. The subject is 
greatly complicated, in fact largely determined, by the fact that entanglement is a very 
fragile resource, in the sense that it inevitably deteriorates while being manipulated or 
stored. We shall see that the attempt to overcome this difficultly has led naturally to 
the consideration of network theory. 

The most direct way to produce entanglement between spatially separated parties 
is to entangle two particles locally and then to send one of them physically to another 
location. Most research in this direction involves sending photons through optical fibres, 
which suffer from inherent limitations due to photon loss via absorption as well as 
decoherence of the state. The limit at this time is about 100 km [28], because the 
probability of transmission decays exponentially with the distance, becoming as low 
as 10 -20 for 1000 km [29J. Following the techniques developed in classical systems 
in order to overcome similar but much less severe limitations, creating entanglement 
between distant stations via a series of intermediate nodes has been proposed. The 
main and crucial difference with classical information is that qubits cannot be copied, 
so that the intermediate links have to be joined together in a very subtle way, known as 
entanglement swapping, which will be discussed below. This one-dimensional system of 
links and nodes has then a natural generalization to a network of arbitrary geometry. 

2.2.1. Local operations and classical communication In manipulating entanglement in 
a quantum system, we typically begin with a given distribution of entangled pairs of 
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Figure 2: Entanglement distillation: the entanglement of two links can be concentrated 
by LOCC into one connection, so that E[o) > max{E(p), E(p')}. 

qubits and then apply a series of operations designed to distribute the entanglement in a 
useful way The problems we consider naturally impose a distinction between local and 
distributed resources. It is important to have a clear notion of what kind of operations 
are allowed on these resources. In this review, we shall only allow local operations and 
classical communication (LOCC) on the network. It turns out that it is possible to 
provide an operational definition of entanglement as the quantum resource that does 
not increase under LOCC. It is easy to see that this definition is equivalent to the pure 
mathematical one given by (12.41) . 

We suppose that initially some entangled pairs of qubits are present in a quantum 
network, as shown in figure dj with one party of each pair occupying a node of the 
network. However, several qubits from different pairs will occupy a single node along 
with other possible resources, both quantum and classical. When we speak of local 
operations and resources, we mean quantum operations and resources within each node. 
On the other hand, we allow only classical messages to be sent between the nodes. 
Qubits belonging to distinct nodes cannot interact quantum mechanically, so that no 
further entanglement can be created between remote stations. However, qubits within 
a node may interact in any way, including with ancillary (local) resources, and any 
measurement may be performed on the components of the node. 

2.2.2. Distillation of weakly entangled states We now address the task of concentrating 
or distilling the entanglement of two (or more) weakly entangled states into a pair with 
higher entanglement, using LOCC only. Suppose that these states are arranged in 
a parallel fashion so that local operations may be performed jointly on all left-hand 
members and on all right-hand members as shown in figure [2J An important task is to 
find criteria determining when and how a given set of states can be transformed into a 
more highly-entangled target state. 

Majorization theory Majorization theory was developed to answer the question: What 
does it mean to say one probability distribution is more disordered than another? One 
application to quantum mechanics is via the connection between disorder (of the Schmidt 
coefficients) and entanglement. Consider an initial pure state \a) and a target state \/3) 
in a bipartite system. These states may have any dimension d, that is, they are pairs 
of qudits. Denoting by cx the unit vector of the d Schmidt coefficients of |a) sorted 
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in decreasing order (and similarly for |/3)), Nielsen showed in [30] that a deterministic 
LOCC transformation from \a) to |/3) is possible if and only if the inequalities 

n n 

J>;<J> (2.13) 

i=0 i=0 

hold for all n G {0, . . . , d — 1}. In this case, cx is said to be majorized by /3, which 
is denoted by cx -< f3; see [31J for more details on this topic. Note, for instance, that 
the maximally entangled state can be deterministically transformed in any other pure 
state. As another example, consider setting p > p' for the two pure states \p) and \p') 
depicted in figure [2J so that one finds cx = (poPo; PoPii PiPoj PiPi)- In this case, the two 
pairs of qubits can be transformed into a single connection \a) if and only if p p' < a . 

Moving now to non-deterministic transformations, the optimal probability for a 
successful LOCC conversion is |25|: 



prob( \a) i y |/3)) = min <J } • (2-14) 

Using this formula, it is trivial to check that a two-qubit pure state \ip) can be 
transformed into a Bell pair with optimal probability 2<p 1; as stated in section 12.1.21 
Explicitly, this result is obtained by performing on one of the qubits a generalized 
measurement defined by the operators 

.^i =a/^|0)(0| 4 |1)(1 




(2.15) 

M 2 = /1-^|0)(0|, 
which is known as the "Procrustean method" of entanglement concentration |32j . 

Distillation The distillation of mixed states is a somewhat more difficult problem 
than that of pure states. Many techniques have been developed for doing mixed-state 
entanglement distillation in connection to quantum error-correction [33J. This work has 
stimulated much research, which is attested by the recent review article [31]. For our 
purpose, however, it is sufficient to notice that: 

(i) In contrast with the case of pure states, at least two copies of a Werner state are 
needed to get, by LOCC and with finite probability, a state of higher fidelity [35] . 

(ii) Perfect Bell pairs can be distilled from N Werner states only in the limit N — > 
oo 1361. 



The first distillation scheme, which was proposed by Bennett et al in [37], is depicted 
in figure HJ the two states p = pw( x ) an d p' = Pw( x ') are distilled into the state 
c = Pw(x"), with 

*" = X + + W (2.16) 
3 + 3xx' K ' 
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Figure 3: Entanglement swapping: the middle station performs a measurement in the 
Bell basis on its qubits. This creates a quantum connection between the two opposite 
stations, which depends on the outcome m of the measurement. 

The resulting state is closer to the target state |$ + ) if both p and p' are entangled (that 
is, if x,x' > 1/3) and if x > x' > 2x/(l + 4x — 3x 2 ). It is important to note that this 
operation is not deterministic since it succeeds only with probability 



This procedure can be iterated, with x increasing after each purification, until it is 
arbitrarily close to 1 (considering perfect operations). The state is thus said to be 
distillable. However, one is often interested in the yield, which is defined as the 
asymptotic ratio of the number of input states to the number of output states. The above 
protocol requires a diverging number of states to produce one arbitrarily pure state, and 
thus has a vanishing yield. In order to get a positive yield, the previous method is applied 
until states of sufficiently large x are generated, and then one switches to distillation 
techniques using one-way communication. Many improvements and variants over this 
construction exist; see [38] for example. 

2.2.3. Entanglement swapping The basic operation to propagate the entanglement in 
a quantum network is the so-called entanglement swapping [TBI EH SOI HI] depicted in 
figure [3j By performing a Bell measurement on the qubits b and b' at the station B, one 
creates a quantum link between the previously unconnected stations A and C. Then, 
a local unitary that depends on the outcome m of the measurement is applied on the 
qubit c, so that the resulting entangled pair has the standard form given in (12.51) or in 
f)2.8p . This operation is equivalent to teleporting b to c. 

For two mixed states p = pw( x ) an d p' = pw(x'), it is easy to see that the 
entanglement swapping produces the quantum state a m = pw(xx') for all measurement 
outcomes, that is, 



In the case of pure states \p) and |p'), however, the result depends on the outcome, 
and one gets either a singlet or a state that is weakly entangled. Remarkably, the average 






V m. 



(2.18) 
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entanglement of the resulting states \a m ) is not less than that of the initial states: 

1 4 

E ^m) = 1*52 E (°m) = 2min{p 1 ,p / 1 }, (2.19) 

m=l 

which is, for |p) = \p'), the "conserved entanglement" described in |42j. This property of 
pure states will be used in the various protocols of entanglement propagation described 
in section O 

Maximizing the average entanglement of the outcomes is of prime importance 
for random or statistical processes, but it is sometimes advantageous to get as much 
entanglement as possible for all outcomes. To this end, we use the rotated Bell basis 
Bx = (X®id2) B to perform the entanglement swapping; see section 1.1 in [43J. In this 
case, one gets four outcomes \a m ) satisfying 



E(a m ) = l-y/l- WpopM Vm, (2.20a) 
or, equivalently, 

C(a m ) = C(p)C(p') Vm, (2.206) 

The entanglement swapping procedure can be iterated, creating a quantum 
connection between qubits that are more and more distant. However, the resulting 
long-distance entanglement that is generated in this way decreases exponentially with 
the number of swappings. This is clear in the case of mixed states and for pure-state 
Bell measurements in the Bx basis; the general proof can be found in section 1.3 in |43j. 
Because of this important loss of entanglement, new schemes have to be designed to 
efficiently entangle any two stations of a quantum network. The various protocols 
proposed so far that achieve this task are described in the following chapters. 

Noisy operations Thus far, we have assumed that all quantum operations are ideal or 
perfect. However, in practice, every operation or measurement introduces an error. Here 
we consider two sources of error: those arising from applying an imperfect gate, (unitary 
operation), and those arising from an imperfect measurement. In each case, we model 
the error by including a small depolarizing term to a channel, which replaces a fraction 
of the density matrix with a completely decoherent (i.e. unknown) state [3,112]. For 
instance, a two-qubit gate Off, with an error probability 1 — 77, acting on a multi-qubit 



state p, is replaced by the map 

p OM = r, Of al [p] + l^L idy ® tr y [p], (2.21) 

where tr^- denotes the partial trace over qubits i and j. In the case of a one-qubit 
rotation, it can be shown that this model is exactly correct for isotropic errors [B|. On 
the other hand, suppose we model a measurement error by assuming that we have a 
small fixed probability 1 — rj of reading 1 when the particle was actually measured into 
the state. In this case the single qubit measurement projectors with errors, 

Po' = v\0)(0\ +(1-77)|1)(1|, (2.22) 

P? = V\W\ + (1 - V) \0)(0\ , (2.23) 
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Figure 4: The nested purification used by the quantum repeaters, (a) Elementary 
connections are continuously generated between neighbouring stations, and an 
entanglement swapping is performed at every second node. The resulting states have 
less entanglement, but they are repeatedly purified by distillation, (b) Once the states 
have regained a sufficiently large fidelity, the procedure is iterated at a higher level. This 
eventually leads to an entangled pair of qubits between the two endpoints of the chain. 



are equivalent to applying the depolarizing channel 

Of[p] = (1 - 2r 1 )p + 77 id! <g> tn[p], (2.24) 

followed by a perfect measurement. Propagating these errors according to this simple 
model allows us to estimate the error resulting from a particular protocol. 

2.3. The quantum repeaters 

A great deal of theoretical and experimental effort has been put forth to distribute 
entangled states over long distances using essentially one-dimensional lattices. As we 
mentioned above, networks were introduced to solve problems caused by unavoidable loss 
and decoherence through free-space and fibre links. In particular, quantum repeaters, 
which have entanglement swapping at their heart, have received the most attention. 

The initial proposals for quantum repeaters use a hierarchical scheme of swapping 
and purification steps jHJ [5j [6] . Including purification in the protocol is necessary once 
one introduces real- world noise and errors. Noise enters the system in two ways. First, 
each operation or measurement will reduce the fidelity of the desired state. This noise 
is modelled by (I2.2ip and f)2.24p . Second, even in the case of perfect operations, if 
one begins with a slightly impure state, then the state resulting from a protocol decays 
rapidly in the number of operations, such as the swapping in (I2.18p . to a useless separable 
state. As explained above, the losses suffered by a state transmitted through the links of 
the network increases exponentially with distance. This introduces a maximum length of 
elementary links because there is a minimum fidelity -F min ( I2.12p below which purification 
is impossible. The repeater protocol first prepares several states of F > -F m ; n along a 
relatively short link and stores them in quantum memories. They are then purified to 
produce a pair with higher fidelity. Entanglement swapping is then performed on two of 
these purified states on neighbouring links thus creating entanglement across a distance 
that is twice the length of the elementary link. The noisy swapping will again reduce 
the fidelity, so that more quantum memories and more purifications are needed. This 
procedure is iterated so that, in principle, highly entangled states can be created across 
long distances. 
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Based on the initial protocol for quantum repeaters, many improvements have been 
suggested; see for a review on this topic. For instance, it has been shown that the 
number of qubits per station does not have to grow with the distance [7]. Yet, the 
realisation of quantum repeaters is still a very challenging task, which is mainly due to 
the need of reliable quantum memories [45] . 

2.3.1. Implementation While many aspects of theoretical quantum information await 
improvements in technology, experiments sending entangled pairs were performed early 
and have seen much improvement. Teleportation, which is the crucial ingredient in 
swapping, has been performed with photons [39, 46J, with recent improvements in the 
direction of telecommunication |47] and atomic systems [48J. 

Entangled photon pairs have recently been separated through air over « 
150 km [49]. While the detected pairs in this experiment are found to be highly entangled 
by violating the CSHS inequality ]50j, it is currently subject to high attenuation rates, 
due mostly to turbulence. One of the first implementations of a quantum network, 
the DARPA Quantum Network, consists of several nodes and supports both fibre-optic 
and free-space links [51]. It is capable of distributing quantum keys between University 
campuses. Another example is the SECOCQ network in Vienna, which successfully 
relayed quantum keys across a network connecting several business campuses [52] . 

Another kind of entangled state providing raw material for entanglement 
distribution schemes can be created as the steady state of a driven, dissipative, 
system ]53j. Proposals have been made to couple a linear system of spins to photons in 
an optical fibre via opto- mechanical transducers [SUES]. By tuning coupling constants, 
a variety of bi- or multi-partite pure entangled states are steady-state solutions to the 
equations describing these systems. But, to date, we know of no work integrating these 
states into entanglement distribution protocols. 
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3. Entanglement percolation 

An approach to entangling distant parties that is conceptually different from the 
quantum repeater was proposed in In this paper, the main question is: given 

a quantum network, which operations should be performed on the nodes so that 
entanglement is best propagated? 

A first answer is presented in section 13.11 for some graphs, the gain obtained from 
distillation of weakly entangled states can balance, or even surpass, the loss of 
entanglement resulting from the swappings. This result is promising, but it is merely an 
adaptation of the quantum repeater protocol to specific graphs and quantum states. In 
fact, it somewhat replaces the repeated generation of elementary links by a deterministic 
accumulation of the entanglement of existing links, provided that they are in a pure 
state, or that the graph has a specific structure. 

The solution proposed in [SB] is of a different nature. Its underlying idea is that in one- 
dimensional networks, any defective link destroys the whole procedure, while in higher- 
dimensional networks the information can still reach its destination through other paths. 
This phenomenon is related to percolation theory, which states that if there are not too 
many defective links in, say, a square lattice, then any two nodes of the network are 
connected by a path with non- vanishing probability. Strategies based on this idea are 
described in section l3~2l Initially limited to pure states, they have since been generalized 
to some special cases of mixed states; see section 13.31 

3.1. Deterministic protocols based on distillation 

In this section, we show that qubits can become entangled over large scales in some two- 
dimensional (planar) graphs, using predetermined sequences of entanglement swappings 
and distillations. Very close to the quantum repeaters in spirit, this method can be 
applied either with mixed states in a restricted class of graphs (section I3.1.ip or with 
pure states in lattices of high connectivity (section f3. 1 .2p . 

3.1.1. Hierarchical graphs Hierarchical graphs iterate certain geometric structures, so 
that at each level of iteration either the number of neighbours or the length of the 
connections increases. Various such graphs were considered in [57]; let us give here yet 
another example in which both pure and mixed state entanglement can be generated 
between nodes that lie at any level of the hierarchy. In this example, two infinite ternary 
trees, in which each link is an entangled state p, are connected at their leaves by a state 
cr; see figure EH The protocol runs as follows: First, two entanglement swappings are 
performed on the central links p, a, and p, yielding a state a. Second, the three states 
a that connect the new leaves of the ternary trees are distilled, leading to a single 
connection a'. The procedure can be iterated as long as E(a') > E(a). In this way, 
nodes lying at higher and higher level of the hierarchical graphs become entangled, that 
is, larger and larger quantum connections are created. In the following, we determine 
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Figure 5: A hierarchical graph that consists of two infinite ternary trees. The links of 
the trees are in an entangled state p, while the central connections may be in a different 
state a. (a) Entanglement swapping is performed at each extremity of the trees, creating 
three links a between every pair of new leaves, (b) Every group of three states a is 
distilled into one connection a'. The process reaches any level of the hierarchical graph 
if E{a') > E(a). 



the minimum amount of pure or mixed-state entanglement of p and a for this strategy 
to be successful. 



Pure states In section 12.2.31 we saw that the result of the entanglement swappings 
depends on the basis that is chosen for performing the two-qubit measurement. In 
order to facilitate the comprehension of the mechanism, we choose the Bx basis, so 
that all outcomes are equally entangled: C(cx) = C 2 (p) C(a); see (12.20 £>[) . The optimal 
distillation of cr® 3 into a', as described in section 12 .2 .2\ requires working with the Schmidt 
coefficient a = — C 2 (a))/2. One then finds a' = max{i, a 3 ,}, and the recursion 

relation for the entanglement of the central connections reads: 

E' = min jl, 2 - i(l + ^1 - /i E(2 - £)) 3 J , (3.1) 

with E' = E(a') and E = E(a). In this equation, the parameter p = C 4 (p) lies in 
the interval [0, 1]. It is an easy calculation to show that E' < E for any value of E if 
p < p c = |, that is, if E{p) < E c « 0.35. This means that the entanglement cannot be 
propagated in the hierarchical graph if the links p are too weakly entangled. In contrast, 
if p > p c , one stable and non-trivial fixed point appears in (j3.ip ; see figure El In this 
case, nodes lying at any level of the hierarchy can be connected by an entangled state 
of two qubits; one further shows that this state is a perfect Bell pair if p > p* w 0.655. 



Mixed states The scenario of a mixed-state hierarchical graph, where the connections 
are Werner states p = pw( x ) an d o = Pw(y)i is quite similar to that of pure states. 
First, two consecutive entanglement swappings are performed on the central states p, 
a, and p, leading to a state a = pw(y) with y = x 2 y. Then, we try to concentrate 
the entanglement of the three states a into one connection a'. As for the pure states, 
we would like this operation to be deterministic, but the distillation of mixed states is 
intrinsically probabilistic. In order to get a result in a predictable fashion, we distill two 
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Figure 6: Graph of the recursion relation (13. ip . Three regimes have to be distinguished: 
/i < fi c (dotted bottom curve), \i c < p, < p* (dashed middle curve), and p > p* (long- 
dashed upper curve). A non-trivial fixed point (bullet) appears only if p > p c , and a 
Bell pair is created after a finite number of iterations if p > p*. 









; 




; 


1 












) 6 6 6 < 












t 


: 




: 


rr. 


: : 


T 



Figure 7: A "centipede" in the square lattice: the entanglement of the links is 
progressively concentrated along the "spine" of the centipede. If the amount of 
entanglement present in the links is larger than a critical value E c , then perfect Bell 
pairs are obtained on the central path and long-distance entanglement can be generated. 



connections only and take the third one if the distillation failed. From ( I2.16P and ( 12.1 7J , 
the distillation succeeds with probability p = (1 + y 2 )/2, and the average entanglement 
of a' = pw(y') is 

2 

y'(x, y) = ^(5 + Ax 2 y - 3x 4 y 2 )- (3.2) 
o 

If the links of the graph satisfy x > x c = a/18/19 and y > y c (x) = (2x — 
Vl9x 2 — 18)/(3a; 3 ), then a stable fixed point y* appears. In this case, iterating the 
entanglement swappings and the distillations generates some long-distance pairs of 
qubits whose entanglement approaches y* = (2x + y/l9x 2 — 18)/(3x 3 ). 

3.1.2. Regular graphs We have just shown that entangled pairs of qubits can be gen- 
erated over a large distance in graphs with a hierarchical structure, under the condition 
that the entanglement of the bonds if larger than a critical value E c . The self-similarity 
of these graphs allows one to design natural sequences of entanglement swappings and 
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distillations but suffers a physical limitation: either the length of the bonds or the num- 
ber of qubits per node grows exponentially with the iteration depth. We now consider 
regular two-dimensional lattices, that is, periodic configurations of links throughout the 
plane in which the nodes have a fixed number Z of neighbours. 

A deterministic strategy to entangle two infinitely distant nodes in lattices in which 
each node has a number of nearest neighbours Z > 4 was proposed in [57]. It is very 
similar to the recursive method developed in the previous example, but it works only 
with pure states. In the case of a square lattice of links \ip), the idea is to sequentially 
shorten the legs of a "centipede", so that the entanglement of the links is gradually 
concentrated; see figure [7J This eventually yields a perfect Bell pair at the spine of 
the centipede, on which infinitely many entanglement swappings can be applied, and 
therefore long-distance entangled pairs of qubits are generated. More precisely, one 
starts by applying two entanglement swappings in the Bx basis on the states \ip), \cr), 
and \(p) at the extremity of each leg. This results in a state \a), and as in the previous 
example, we have C(a) = C 2 (ip) C(er). The difference is that the distillation is now 
performed on \a) <8> \f) rather than on \a)® 3 . Very similarly to (13. ip . one finds the 
following recursion relation: 



where \x = C A ((p) is a function of <pq. One can show that there always exists a non-trivial 
stable fixed point for this equation. However, the fixed point of (13.31) is strictly smaller 
than unity when E((p) < E c rs 0.649. In this case, although we do concentrate some 
entanglement along the spine of the centipede, we still face the problem that the spine 
is a one-dimensional system, which therefore exhibits an exponential decrease of the 
entanglement with its length. On the other hand, if E((p) > E c , then the fixed point is 
reached in a finite number of steps, and a maximally entangled state is generated. Since 
the spine now consists of perfect connections, any two nodes lying on it can share a Bell 
pair, regardless of their distance. 

3.2. Percolation of partially entangled pure states 

We have demonstrated that a way to generate long-distance entanglement in a lattice is 
to distill a "backbone" of Bell pairs and then to perform some entanglement swappings 
along this path. Three conditions have to be satisfied for this method to work. First, 
the nodes that belong to the backbone must have at least four neighbours each: two 
connections are part of the backbone, while the other two are used to distill the former. 
Second, the entanglement of the bonds has to be larger than a critical value E c that 
depends on the lattice geometry. Third, the links have to be pure states and not mixed 
states, because the distillation of a finite number of Werner states never leads to a Bell 
pair (section |2.2.2p . Since this deterministic strategy creates a chain of Bell pairs by 
using only a strip of finite width from the lattice, it seems that it does not exploit 
the full potential of two-dimensional networks. In this section, we review the method 




(3.3) 
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of entanglement percolatioiig] that was published in [56J and that partially relaxes the 
above conditions: 

(i) percolation is a genuine two- (or multi-) dimensional phenomenon, and thus it 
applies to any lattice; 

(ii) entanglement percolation undergoes a phase transition with the entanglement of the 
links, but the corresponding critical value is smaller than that of the distillation 
method [561 E3 [60]; 

(hi) in certain lattices, this method applies also to two-qubit mixed states of rank less 
than four |6T1 

The simplest application of entanglement percolation in infinite lattices is presented 
in section 13.2.11 In this case, the connection to classical percolation theory is 
straightforward and entanglement thresholds are readily determined. Then, we show 
that modified versions of this classical entanglement percolation (CEP) yield lower 
thresholds. For instance, we reconsider the hexagonal lattice with double bonds 
proposed in [56], in which quantum measurements lead to a local reduction of the 
entanglement but change the geometry of the lattice. This operation increases the 
connectivity of the graph and lowers the entanglement threshold, an effect which is called 
quantum entanglement percolation (QEP); see section 13.2.21 Finally, an alternative 
construction that uses multipartite entanglement is given in section 13.2.31 Using 
incomplete measurements, this protocol creates entangled states of more than two qubits 
and improves not only the entanglement threshold, but also the success probability of the 
protocol for any amount of entanglement in the connections for all the cases considered 
in section 13.2.31 



3.2.1. Classical entanglement percolation Classical percolation is perhaps the funda- 
mental example of critical phenomena, since it is a purely statistical one [53] ■ At the 
same time, it is quite universal because it describes a variety of processes, with ap- 
plications in physics, biology, ecology, engineering, etc. [64). Two types of models are 
typical considered, site- and bond percolation. In bond percolation, the neighbouring 
nodes of a lattice are connected by an open bond with probability p, whereas they are 
left unconnected with probability 1 — p; see figure In site percolation, the sites (i.e. 
nodes) rather than bonds are occupied with probability p. In either case, for an in- 
finite lattice, one would like to know whether an infinite open cluster exists, that is, 
whether there is a infinitely long path of connected nodes. It turns out that an unique 
infinite cluster C appears if, and only if, the connection probability p is larger than 
a critical value p c that depends on the lattice. Few lattices have a threshold that is 

f Note that the probabilistic nature of quantum physics makes percolation theory a particularly well- 
adapted toolbox for the study of quantum systems which undergo, for instance, measurements. Ideas of 
percolation theory are for example useful in the context of quantum computing with non-deterministic 
quantum gates [5H]. We refer the interested reader to pp. 287-319 in [59] for an overview of the 
application of percolation methods to the field of quantum information. 
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(a) p = 0.25 (b) p = p c = 0.5 (c) p = 0.75 



Figure 8: In classical bond percolation, the bonds are open (present) with probability 
p and closed otherwise. Groups of nodes connected by open bonds are called clusters. 
In these examples, the largest cluster is highlighted: typically, it is small and bounded 
for p < p c , but it spans a finite proportion of the lattice for p > p c . In the latter case, 
there is a unique infinite cluster C. 



exactly known. Among them, we find the important honeycomb, square, and triangu- 
lar lattices, withpP = 1 — 2 sin(7r/18), p° = 1/2, andp^ = 2sin(7r/18), respectively [65J. 

Suppose now that we want to generate some entanglement between two distant sta- 
tions A and B in a quantum lattice, where each connection denotes a partially-entangled 
pure state \<p). Classical entanglement percolation (CEP) runs as follows [3B]: every 
pair of neighbouring nodes tries to convert its two-qubit state into a Bell pair, which 
succeeds with an optimal probability p = E(ip); see section [2.2.21 If this value is larger 
than the threshold p c of the lattice, that is, if the entanglement of the links is large 
enough, then an infinite cluster C appears. The probability that both A and B belong 
to C is strictly positive, and in this case a path of singlets between these two nodes can 
be found. Then, exactly as described in the previous section, one performs the required 
entanglement swappings along this path such that A and B become entangled. Note 
that the path of Bell pairs is randomly generated by the measurement outcomes at the 
nodes, which contrasts with the deterministic location of the "backbone" generated by 
the distillation method. 

A quantity of primary interest when studying the efficiency of CEP is the percolation 
probability 

0(p) = P{AeC), (3.4) 

which is the probability that a node A belongs to the infinite cluster. This value is closely 
related to the percolation threshold: in fact, in an infinite lattice we have 8(p) =0 for 
p < p ci whereas 8(p) > for p > p c . In our case, we are interested in the probability 
P(Ao-o B) of creating a Bell pair between two nodes A and B separated by a distance L. 
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For p < p c , this probability decays exponentially with L/£(p) [66], where the correlation 
length £(p) describes the typical radius of an open cluster. Above the critical point, the 
two nodes are connected only if they are both in C. In the limit of large L, the events 
{A G C} and {B G C} are independent, so that 

P(AooB) = 6 2 (p), (3.5) 

and therefore the problem is reduced to studying 9{p). 

3.2.2. Quantum entanglement percolation A natural question is whether the 
entanglement thresholds E c = p c defined by the classical percolation theory are optimal. 
In fact, percolation of entanglement represents a related but different theoretical 
problem, where new bounds may be obtained. This is of course equivalent to determining 
if the measurement strategy based on local Bell pair conversions is optimal in the 
asymptotic regime. Several examples that go beyond the classical picture were developed 
in [551 EH ED], proving that CEP is not optimal; such results are referred to as "quantum 
entanglement percolation" strategies. All these examples are based on the average 
conservation of the entanglement after one swapping, as described in section 12.2.31 but 
they do not provide a general construction to surpass the classical percolation strategy. 
In the following paragraph we review the original example of [56], since it gives some 
insights into the way a quantum lattice can be transformed by local measurements, 
and we let the reader consult the articles [571 ED] for the other examples. Finally, an 
improved strategy that makes use of multipartite entanglement will be described in 
section 13.2.31 



Honeycomb lattice with double bonds Let us consider a honeycomb lattice where each 
pair of neighbouring nodes is connected by two copies of the same state \<p); see figure l9al 
The CEP protocol converts all bonds shared by two parties into a single connection, 
and from majorization theory we know that a double bond can be optimally distilled 
into one pair of qubits with entanglement E' = 2(1 — cp*). Setting this value to be equal 
to the percolation threshold for the honeycomb lattice, one finds that the entanglement 
can be propagated if E(ip) > E c , with 

E c = 2 ^1 - yi-^j « 0.358. (3.6) 

However, there exists another measurement pattern yielding a better percolation 
threshold (figure l9b|) : half the nodes perform on their qubits three entanglement 
swappings, which maps the honeycomb lattice into a triangular one. Since the 
entanglement of the connections is not altered, on average, it follows that a lower 
threshold is found: 

E c = p*& 0.347. (3.7) 

This proves that CEP is not optimal since it cannot generate long-distance entanglement 
in the range E((p) G (E C ,E C ), whereas the quantum entanglement percolation achieves 
it with a strictly positive probability. 
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(a) (b) 

Figure 9: Each bond of the honeycomb lattice consists of two copies of the state \<p). 
(a) The dashed nodes perform three entanglement swappings in the Bell basis B. (b) A 
triangular lattice of identical average entanglement is obtained, and CEP is now possible 
in the new lattice. 



3.2.3. Multi-partite entanglement percolation We have seen that CEP can be enhanced 
by first applying some quantum operations at the nodes [56j EH EQ]. All examples 
proposed in these articles consist of transforming the quantum lattices by a sequence 
of entanglement swappings, thus conserving the average entanglement of the bonds. 
They are, however, restricted to purely geometrical transformations, and they apply 
to specific lattices only. In this respect, it was not clear whether the CEP strategy, 
and particularly the corresponding threshold, could be improved in general. A positive 
answer to this question was given in [67], in which a class of percolation strategies 
exploiting multi-partite entanglement was introduced. To that end, one performs the 
entanglement swappings in a more refined way, which we describe below. 



Generalized entanglement swapping The key ingredient of the multi-partite method is 
to consider a generalized entanglement swapping at the nodes: First, an incomplete 
measurement (i.e., not a complete projection) is performed at the central node by 
applying the operators 

M 1= |0)(00| + |1)(11|, 



M 2 = |0)(01| + 1 1X10| 



(3.8) 



for which the completeness relation Y^h=i ^l^i = i°U i s satisfied. This measurement 
leaves a two-dimensional subspace at the central station entangled with the two outer 
nodes. Thus, the central station still plays a role in the propagation of entanglement 
through the lattice. Second, not only two but n > 2 links are "swapped" at the same 
time, and then the Procrustean method of entanglement concentration is performed 
on the resulting state; see ( I2.15p . These operations succeed with a finite probability 
that depends on the number of links and on their entanglement [67], generating the 
Greenberger-Horne-Zeilinger (GHZ) state 

|GHZ n ) ee |U; , (3.9) 
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Figure 10: Multi-Partite entanglement percolation, (a) Some nodes (filled circles) of 
the quantum lattice £ probabilistically transform two links into a GHZ state on three 
nodes, which is depicted by a triangle. No operation is performed at the other nodes 
(empty circles), (b) The nodes of the new lattice £ represent the GHZ states that have 
been created, and a bond is present whenever two GHZ states are adjacent in £. 



which is the generalization of the Bell pair |$ + ) to n qubits. 

Entanglement thresholds: from bond to site percolation In the multi-partite strategy 
one creates from the links \(p) of a quantum lattice £ a new lattice £, where the nodes 
represent the GHZ states obtained with probability p by the generalized entanglement 
swappings. Two vertices in £ are connected by a bond if the corresponding GHZ states 
share a common node in the original lattice. This defines a site percolation process with 
occupation probability p, and above the site percolation threshold of the new lattice, 
the entanglement is propagated over a large distance as follows. Consider the situation 
in which two GHZ states of size n and m sharing one node have been created. One 
builds a larger GHZ state on n + m — 1 particles with unit probability by performing 
a generalized entanglement swapping on the two qubits of the common node. This 
operation is iterated, eventually yielding a giant GHZ state spanning the lattice. Then, 
given a GHZ state of any size, a perfect Bell pair is created between any two of its qubits 
by measuring all other qubits in the X basis. 

An example of multi-partite entanglement percolation is given in figure [TUl In this 
case, since the probability to create a GHZ state of three qubits is equal to that of 
creating a Bell pair (only two links are required), the minimum amount of entanglement 
of the bonds for generating long-distance entanglement is E c pa 0.650, which is equal to 
the site percolation threshold in £ [6B]. Since the critical value for CEP is E c ~ 0.677 
(the threshold for bond percolation in £ [69]), it follows that multi-partite entanglement 
percolation surpasses CEP in the range E(<p) G (E C ,E C ). 

The previous example appeared in [67] together with many other lattices for which 
the critical values using the multi-partite strategy are lower than the thresholds for CEP. 
In particular, an improvement is found for all Archimedean lattices, which are tilings of 
the plane by regular polygons (such as the square or the triangular lattice). Moreover, it 
is shown that not only the thresholds but also the probabilities P(Ao-oB) are better for 
any amount of entanglement E c < E(tp) < 1. This clearly indicates that the interplay 
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of geometrical lattice transformations and multi-partite entanglement manipulations is 
a key ingredient for propagating entanglement in a quantum network. 

3.3. Towards noisy quantum networks 

We have already noted that creating perfectly entangled states via LOCC by consuming 
a finite number of states on a network of links of full-rank mixed states is not possible, 
even considering perfect operations. Two directions we may take from here are: i) If we 
cannot create a singlet, we may try to create a state with the highest possible fidelity, 
ii) Ask instead: For what class of mixed states it indeed possible to create a singlet, and 
what are the optimal protocols? 

We have already paid some attention to the first question above and will review a 
more detailed examination in section 15.41 The answer to the second question is based 
on two results: i) It is not possible to transform a single copy of a mixed state to a 
pure state with local operations, and ii) In the case of two-qubit pairs, a pure state 
may only be obtained from two or more pairs belonging to a certain class of rank-two 
states [70], which we shall call purifyable mixed states (PMS). Protocols on networks 
using this approach were studied in pjTj 16*2] . This work was extended to a hybrid 
approach addressing both the first and second question in [71] . The protocols presented 
in these three articles are the subject of this subsection. Note that what these authors 
call a PMS is a slightly more restricted class of states. 

The most obvious constraint on the design of entanglement distribution protocols 
is that at least two disjoint paths of PMSs must exist between two stations in order to 
have a non-zero probability of a singlet between them. The final stage must consist of 
purifying multiple PMSs. We optimally convert two PMSs to a singlet in two stages. 
First, we perform a pure-state conversion measurement (PCM) as follows. We perform 
a quantum logic operation consisting of unitaries called the controlled-NOT gate at each 
local station, with qubits from one pair acting as targets in both cases. We then measure 
the targets in the computational basis, and if both results are 1, we have generated a 
pure state. If the two input states were identical, then the output state on success is 
already a singlet. But, in general, another singlet conversion using the Procrustean 
method according to (12.151) in section 12.2.21 must be performed. However, as we shall 
see below, it is sometimes advantageous to delay this final singlet conversion and use 
the intermediate state in a different way. 

Swapping The behaviour of PMSs under swapping is similar to that of pure states. 
As in the case of pure states, we project onto the Bell basis, but now only two of the 
resulting states are useful, themselves being PMSs. Most importantly, the fidelity of the 
average resulting state decays exponentially in the number of links as it does for pure 
states. 
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CEP CEP protocols analogous to pure state protocols are denned by taking regular 
lattices with multiple PMSs per bond. This is the most basic way to provide the 
necessary two pairs between nodes. This situation is similar to the double-bond 
honeycomb lattice of section I3.2.2[ except that the states are PMSs and there may be 
more than two pairs connecting nodes. We saw for pure states that a singlet conversion 
on each bond succeeds with a probability p, mapping the entanglement distribution 
problem to classical percolation with bond density p. In the present case, we still map 
directly to classical percolation, but the bond density p is determined by the success 
rate of some conversion protocol of the PMSs to a singlet. For two PMSs, the optimal 
protocol for this conversion is known [70] and has a maximum success rate p = 1/2, 
so that, for instance, percolation is possible on the double-bond triangular lattice but 
not on the double-bond square lattice. It is possible to achieve p > 1/2 for three 
or more bonds, but the optimal conversion protocol is not known in this case. One 
protocol projects locally the entire state onto the subspace that is pure and entangled— 
the distillable subspace [72]. This was investigated in [62], along with better protocols 
that distill multiple bonds in smaller groups and reuse states from some of the failed 
conversions. 

QEP We have seen that for pure-state entanglement percolation, the first step in 
quantum pre-processing that goes beyond CEP is to note that, for a chain of two pairs, 
swapping before singlet conversion is better than singlet conversion before swapping. 
For PMSs the simplest analogy has more choices of when to do a pure-state conversion 
(PCM), singlet conversion, or swap, because we need a minimum of two pairs in parallel 
for each link in the chain. It turns out that the optimal method for this case of two pairs 
per link is to perform the PCM on each link and then to swap the resulting pure states 
before doing a singlet conversion. If the input states are identical, this is equivalent to 
CEP because, as mentioned already, a successful PCM already returns singlets. In [62] 
this method was applied to small configurations and the results, in turn, to some regular 
and hierarchical lattices. 

In all of these protocols, the only non-trivial case, with respect to CEP, is when the 
multiple pairs making up a link are PMSs with differing parameters; otherwise we get 
CEP again after the first purification step. In the case where the parameters are different 
within a link, we must first do a PMS, which maps the problem to the original pure-state 
percolation problem with some of the links probabilistically deleted. A similar idea is 
used in |71J, where pairs in certain rank-three states replace the PMSs. These can be 
converted to binary states via a sort of PCM that leave a separable state on failure. 
The resulting lattice can be treated with error correction methods as in section HJ with 
the difference that the pairs in binary states are deleted probabilistically. 
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3.4- Open problems 

It is natural to wonder about the optimality of the protocols based on entanglement 
percolation, either for pure states or for certain classes of mixed states. We focus on the 
square lattice here because it is very common, but the arguments apply to other lattices 
equally well. 

It is obvious that percolation protocols cannot be optimal for every amount of 
entanglement of the links. In fact, we have seen in section 13.1.21 that a long-distance 
perfect Bell pair can be obtained deterministically if E(tp) is larger than the threshold 
E c rs 0.649, whereas this is possible in entanglement percolation only if E(ip) = 1. On 
the other hand, the distillation method completely fails if E(<p) < E c , while CEP yields 
positive results in the range E((p) G (|, E c ). Consequently, considering one strategy only 
for all values of entanglement is not sufficient, but finding the optimum one for a given 
value E(ip) in the links is a formidably difficult problem \57\ IB"0"t 157] . In this respect, 
multi-partite entanglement percolation is well-suited to generate long-distance quantum 
correlations regardless of the entanglement of the links, since it leads to high connection 
probabilities and low thresholds at the same time. A somewhat more tractable question 
is: 

Does there exist a value of entanglement per link below which it is impossible to 
entangle two infinitely distant qubits using LOCC in a two-dimensional quantum 
network? 

In the case of pure states the answer to this question is not known [331 EZj. Strategies 
based on multi-partite entanglement percolation are among the most efficient ones that 
achieve this task, but other efficient protocols may exist. For instance, percolation, 
distillation and swapping have been mixed in a hybrid protocol that yields an 
improvement in settings that were resistant to attempts employing distillation and 
swapping separately. This hybrid protocol works as follows. The first step in the 
entanglement swapping protocol is not altered, that is, a Bell measurement on the 
central qubits is performed. There are essentially two entanglement outcomes for the 
two remaining qubits: either a Bell pair or a partially entangled pair is produced. In 
the usual swapping procedure, the process terminates if the former state is produced 
or, in the event of the latter state, a singlet conversion is performed. In the hybrid 
procedure, a partially entangled outcome is instead used together with an existing 
parallel bond in a distillation. For example, using this protocol, a triangular lattice 
is converted to a honeycomb lattice with enough entanglement per link to overcome 
the geometrical disadvantage. Another example involves the hybrid protocol and the 
generalized entanglement swapping applied to the square lattice, which results in the 
lowest known percolation threshold for quantum percolation on this lattice; see chapter 
2.4 in [33] • Finally, better protocols using the quantum network in a completely different 
manner may exist. For instance, it was shown in [73] that long-distance quantum 
correlations can be obtained in a square lattice using techniques of error correction 
(section BSD) - 



Entanglement percolation 



30 



Quite surprisingly, the situation turns out to be opposite in the case of mixed-state 
networks. In fact, suppose that the connections of the network are given by the Werner 
state pw(p) — P \^ + ){^ + \ + (1 ~p) icU/4, with p smaller than the (classical) threshold p c 
for bond percolation in the corresponding lattice. That is, the quantum state describing 
the whole system is a classical mixture of lattices whose links are either perfect Bell pairs 
or completely separable states. In the limit of infinite size, however, none of these lattices 
possesses an infinite cluster of Bell pairs. The threshold p c is thus a lower bound on p 
for a lattice of states pwip) since, by definition, no local quantum operation can create 
entanglement from separable states. In the square lattice, for instance, genuine quantum 
correlations cannot be generated over arbitrarily large distances if p < p° = 1/2, even 
though all connections are entangled in the range p E (|, |). Finally, the situation for 
mixed-state lattices in three dimensions is similar to that of pure states: there exists a 
threshold to generate long-distance entanglement in the cubic lattice (section l4.2.3p . but 
no interesting lower bound is known. In fact, the previous argument leads to a lower 
bound given by the percolation threshold on the cubic lattice p c ~ 0.2488 [74J, but at 
this bound the quantum connections are useless in any case since pwip) is separable 
for p < |. The previous argument can be generalized to any mixed state using the 
concept of the best separable approximation (BSA) to an entangled state, introduced 
in [75]. Given a quantum state p, one decomposes it as the mixture of an entangled and 
a separable state, pE and ps, with positive weights ps and pe such that ps + Pe = 1: 

P = Ps Ps+Pe Pe- (3.10) 

The BSA to the state p is defined by the decomposition maximizing p$. Clearly, if the 
states in a network are such that the separable weight of its BSA ps is smaller than the 
network percolation threshold, there is no protocol allowing long-distance entanglement 
distribution. 

In conclusion, it is of great interest to determine if there exists a lower bound E min 
for propagating the entanglement in quantum networks, and, if this is the case, to design 
new protocols that bring E c as close to E min as possible. 
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4. Network-based error-correction 

In the previous section, we showed that percolation allows one to efficiently create 
entangled pairs of qubits over a large distance in quantum networks that consist of 
pure states or of a restricted class of mixed states. This is a considerable improvement 
over one- dimensional systems, in which the probability to generate remote entanglement 
decreases exponentially with the distance. The connectivity of the nodes plays a key 
role in this respect, but it is not clear if a similar effect exists in networks subject to 
general noise, that is, described by mixed states of full-rank. In fact, a finite number of 
such states cannot be distilled into maximally entangled states (section l2.2.2p . which is 
an essential requirement for entanglement percolation. 

The aim of this section is to review the strategies that have been developed for 
propagating the entanglement in quantum networks whose connections are full-rank 
mixed states. However, we restrict our attention to Werner states defined in ( 12. 8p . This 
entails little loss of generality because any mixed state can be transformed to a Werner 
state with the same fidelity via depolarization. 

4-1. Correction of local errors from a global syndrome pattern 

At present, few schemes have been proposed to generate long-distance quantum 
correlations in noisy networks [?3j [?6j EH UM Although the quantum operations 
that are performed at the nodes are quite different for each scheme, the underlying 
principle is similar: 

(i) The bonds are used to create a multi-partite entangled state that is shared by all 
nodes of the network. Due to the noise in the system, the generation of this state 
is imperfect; 

(ii) Local measurements on all but two distant qubits partially reveal at which places 
the noise altered the creation of the multi-partite state; 

(iii) A global analysis of the measurement outcomes determines the operations that 
have to be applied on the remaining two qubits in order to get useful remote 
entanglement. 

These steps are described in more detail in what follows, and the schemes are reviewed 
in section 14.21 

4-1.1. Creation of a multi-partite entangled state A first hint of the usefulness of multi- 
partite quantum states was given in section 13.2.31 In that setting, a giant GHZ state 
is created on the lattice by extracting perfect entanglement from the bonds adjoining 
each node. Then, measurements in the X basis of all but two qubits transform the giant 
state into a Bell pair between the two remaining qubits. Finally, a local basis rotation 
depending on the measurement outcomes further converts the Bell pair into, say, the 
maximally entangled state |$ + ). The procedure is rather similar here, but the links 
of the networks are used to create either a GHZ state (section 14.2. ip . a surface code 
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(section T4.2.2p . or a cluster state (section I4.2.3p . The key point of the construction is 
that while these states are simple enough to be created by local operations on the nodes 
of the quantum network, they also are tolerant of a certain amount noise, as described 
in the following. 

4-1-2. Syndrome pattern In the protocols involving pure states, it is known exactly 
whether a conversion of partially entangled states into Bell pairs succeeds or fails, since 
this information is given by the outcome of a measurement. In contrast, in mixed-state 
quantum networks, some noise enters the system randomly, and there is no way, a priori, 
to know where this happens. In fact, because every connection is a Werner state with 
non-unit fidelity, the generation of the multi-partite state based on such connections 
leads to a quantum state that contains errors (for instance, bit-flip and phase errors on 
some of its qubits). Hence, if one decides to measure all but the two target qubits right 
after the generation of the multi-partite state, then the choice of the final basis rotation 
will be correct with a probability of approximately only one fourth. This means that 
we have no knowledge at all about which one of the four Bell pairs we are dealing with, 
or in other words, the qubits are in a separable quantum state. 

One strength of the network-based error correction is that one can gain some 
information about the errors without damaging the long-distance quantum correlations. 
In fact, the multi-partite entangled states created in the quantum network are highly 
symmetrical and satisfy a set of eigenvalue equations that can be checked by local 
measurements: if the outcomes do not match the symmetry of the target state at the 
corresponding nodes, which is called a syndrome, then one immediately knows that at 
least one adjacent link inserted an error into the system; see section fl~2l The question 
of determining which link is responsible for the syndrome is treated in what follows. 

4-1-3. Error recovery Syndromes are defined on the nodes of the dual lattice of the 
quantum network, which is either a square lattice (first two schemes) or a cubic lattice 
(third scheme). The generation of the multi-partite entangled state is such that the 
noise entering the system corrupts every link of the dual lattice with probability p. 
This creates chains of errors, which are consecutive corrupted links of the dual lattice. 
Syndromes correspond to the endpoints of these chains and thus come in pairs, as 
depicted in figure lllal 

If one knew the location of all chains of errors, then it would be possible to perfectly 
restore the target multi-partite state. The difficulty of the error recovery is that while 
the positions of the syndromes are known, no other information about the chains is 
available. Since different chains of errors can lead to a similar syndrome pattern, the 
recovery is ambiguous and may lead to a wrong correction of the errors. However, 
Dennis et al. showed that, in an infinite square lattice, a (partial) recovery is possible 
if the error rate p does not exceed a critical value p c [80J. This threshold is found via a 
mapping to the random-bond Ising model and is approximately equal to 10.94%, which 
is numerically calculated in [8*T| . In order to obtain long-distance quantum correlations 
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Figure 11: Error correction from a pattern of syndromes, (a) As the multi-partite 
entangled state is created from the noisy bonds of the quantum network (not shown), 
chains of errors (red lines) appear in the dual lattice (here, a square lattice), and 
syndromes (filled circles) are defined as the endpoints of the chains, (b) In dashed 
blue lines, a possible minimum-weight perfect matching of the syndromes. Inferring 
chains of errors that do not match the actual ones is equivalent to adding defective 
links, and closed loops of chains of errors are created, (c) The error correction fails if 
closed loops extend from one border to the opposite one. In an infinite square lattice, 
this happens if the error rate p is larger than the critical value p c ~ 0.11. 



for p < p c , one should be able to compute all patterns of errors that lead to the measured 
syndromes and then to choose the one that most likely occurred. This is in-feasible in 
practice, but for small error rates p, a good approximation of the optimal solution is 
the pattern in which the total number of errors is a minimum. In fact, such a pattern 
may be efficiently found by a classical algorithm, known as the minimum-weight perfect 
matching algorithm jS2J [S3] . Illustrations of these ideas are given in figure lllbl and 
figure lllcl 

4-2. Examples of protocols 

Now that the general concepts about network-based error correction have been 
introduced, we describe in the following paragraphs the quantum operations that are 
required to generate entanglement over long distance in noisy networks, as proposed 
in [751 US EZ] ■ 

4-2.1. Independent bit-flip and phase errors Any entangled state of two qubits can 
be transformed by LOCC to the Werner state defined in (12.8p . but to understand the 
scheme of [73] it is more appropriate to consider a slightly different parameterization of 
a two-qubit entangled mixed state: 

p(e b ,e p ) = ((1 - e b )(l - e p ), e b (l - e p ), e p (l - e b ), £b£^j ■ (4.1) 

In this equation, e b and e p stand for the probability that the second qubit of the ideal 
connection |$ + ) has been affected by a bit-flip and a phase error, respectively. This 
state is as general as a Werner state in the sense that any quantum state of two qubits 
can be brought to this form using LOCC only. 
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Figure 12: Generation of a giant GHZ state. Every station is connected to its neighbours 
by four entangled pairs of qubits. These links are consumed, via generalized Bell 
measurements followed by local rotations, to generate a GHZ state that spans the whole 
network. An auxiliary qubit (c) is used to check the parity of the "incoming" qubits a 
and b: after two controlled-NOT gatedjjt is measured in the computational basis. The 
station is then tagged with the outcome or 1 of the measurement, which will be used 
to create the pattern of syndromes. 



Protocol in the case of bit-flip errors only The links of a iV x N square lattice are used 
to create a giant GHZ state of N 2 qubits; see figure [T2J For each direction, the stations 
perform the generalized entanglement swapping described in (13. 8p . If one temporarily 
assumes that phase errors are not present in the links of the network, then the resulting 
multi-partite state is a mixture of GHZ states whose qubits are flipped with probability 
e 6 . If the error rate does not exceed the critical value p c ~ 10.9%, then most bit-flip 
errors can be corrected, as depicted in figure [T3"l 

The bit-flip error correction presented above can be applied to arbitrary planar 
networks, as shown by Broadfoot et al in [71]. They also prove that it can be generalized 
to entangled mixed states of rank three, but another method has to be used in the case 
of full rank mixed states, which is the topic of the following paragraph. 



Protocol including both bit-flip and phase errors The global error correction works 
exactly as described above, but each qubit is replaced by a logical qubit which is an 
encoded block of n qubits. Furthermore, all quantum operations on the logical qubit 
are implemented by an appropriate protocol at the encoded level. Phase errors are then 
suppressed by the redundancy of the following code: 

|o>-Ho) = 4f(|+)°*+ l-> 0n )> 

f (4-2) 

ii)^io) = — i-r), 

f Under a controlled-NOT operation, the state 1 11) becomes |10), 1 10) becomes 1 11), while 1 00) and 
1 01} are unchanged. 
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Figure 13: Syndrome pattern for the giant GHZ state, (a) Due to the sequence of 
measurements detailed in figure [12], links (thick red lines) that insert a bit-flip error into 
the GHZ state invert the value of all rightward and upward parity checks, (b) Syndrome 
(circles) are created when an unit cell is surrounded by an odd number of parity check 
"l"'s. The minimum weight perfect matching of these syndromes in the dual lattice 
(dashed lines) reveals the location of the noisy links. 

where |±) are the eigenvectors of the X basis. Using majority vote, up to t phase errors 
can be corrected for each block of n = 2t + 1 qubits. A detailed discussion of a fault- 
tolerant implementation of this encoding is done in [73] , but here it is sufficient to state 
the final result: long-distance quantum correlations can be obtained with the encoded 
version of the global bit-flip correction, but the number of links between neighbouring 
stations has to increase logarithmically with the distance. The main difference with the 
one-dimensional quantum repeaters is that all operations can be applied simultaneously, 
so that no quantum memory is needed. 

4-2.2. Surface codes Surface codes have been used to implement error correction in 
quantum computation and communication jSU EH ESI E7J ESI EH] ■ Here, we review three 
ways in which the surface code is used to generate entanglement over large distances. 
The protocol described in [75] differs from the previous protocols in that the data to be 
transmitted is encoded across the network as the multi-partite state is created. In [78J, 
the authors propose to use the entanglement existing in the links of the network to 
apply stabilizer measurements and to create the surface code across the network. In 
both proposals the number of required entangled pairs scales with the distance at which 
the entanglement is to be generated. In what follows we describe these two protocols, 
whereas the third one, which requires a constant number of quantum connections only, 
is described in section 14.2.31 

Planar surface code For the planar surface code, half of the qubits encode the data and 
half are only used to measure syndromes. The code is implemented via stabilizers. A 
stabilizer of \ip) is an operator A satisfying A = Data is encoded by preparing 
the array in a simultaneous eigenstate of a set of commuting stabilizers, which are defined 
as follows. A lattice is formed by associating an edge with each data qubit. Each face 
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Figure 14: Surface code of depth d encoding a logical qubit. Data qubits are open circles. 
Syndrome qubits are filled circles. To guide the eye, four of the stabilizers are shaded 
grey. Two of the shaded stabilizers are on a boundary, so they have only three Pauli 
operators. The six steps in measuring a stabilizer are shown in detail: A preparation 
measurement, four CNOT gates, and final measurement. Adapted from A.G. Fowler et. 
aL 2010. 



represents a stabilizer formed by the tensor product of four Pauli Z operators tensored 
with the identity on all other data qubits. Vertex stabilizers are defined similarly, 
substituting the X operator for the Z operator. The exceptions are on the edges of the 
array, where the stabilizers are products of only three Pauli operators. The dimensions 
of the lattice are chosen to give two boundaries with partial vertex stabilizers, and two 
with partial face stabilizers, which ensures that the dimension of the space satisfying 
the stabilizer conditions is 2. That is, the array encodes one logical qubit. Details of 
surface codes and error correction can be found in [SU] . 

Surface code communication protocol The planar surface code is implemented on a two- 
dimensional rectangular array of entangled qubits and encodes a single logical qubit; 
see figure [EH We begin with a long rectangular array that supports this code. Square 
sections on the left and right ends also each support a surface code. The left end is 
prepared in a logical state \^/) and the right end in the logical state |0). A sequence of 
operations is performed in such a manner that it spreads the state until the entire 
array encodes \4f). Finally, we perform some measurements that leave the array on the 
right end in the state |^). Errors are measured during the entire procedure and passed 
classically to the right end, where they are processed to correct the result. 

Implementing the surface code communication protocol The array is initialized with 
on the left and |0) on the right, while all the data qubits in the middle are measured in 
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Figure 15: Repeater based surface code. The middle section of the array of qubits is 
split at each syndrome qubit by coupling to a Bell pair. Adapted from A.G. Fowler et. 
al, 2010. 



the Z basis. The state \^f) is then spread across the array by measuring each stabilizer a 
number of times equal to code depth d, which is the same as the height of the lattice. In 
general, increasing d allows better correction for gate errors. The sequence of six gates 
shown in figure [Til is designed so that all stabilizers may be measured simultaneously. As 
a result of the repeated measurements, the syndromes mark chains of errors that extend 
in time (for d discrete steps) as well as space. A change in a stabilizer measurement 
signals a syndrome marking the beginning or end of a chain. These error chains can be 
corrected in the same way as described in the previous subsection. 

As it stands, this protocol creates and transports entanglement only through local 
interactions. Thus, it is not sufficient for communication between distant parties. For 
long distance communication, each syndrome qubit is replace by two syndrome qubits, 
each of which is coupled to one party of a Bell pair as shown in figure [151 The Bell pairs 
must be generated for each of the d steps in spreading the surface code. 

Assuming no loss in transmission and that all gates within the repeater nodes are 
perfect, the average time to failure of a single link only grows non negligibly with the 
code depth d if the fidelity of the entangled pairs satisfies F > 0.92. Furthermore, 
the authors find that by modelling loss in transmission as measurement in an unknown 
basis, loss rates less than 0.45 can be handled efficiently They find that realistic gate 
error rates do not affect these results significantly. The number of qubits per repeater 
grows with the code depth d, which scales as log(N/p c ), where N is the number of links 
and p c is the desired communication error rate. Reference [76] includes a more thorough 
discussion of the error rate and transmission rate. 

Network-based surface code protocol The authors of [78] suggest a different scheme 
to generate long-distance entanglement based on the surface code. The main idea 
is to use the entanglement present in the links of the network to perform stabilizer 
measurements on qubits sitting in different nodes. More specifically, the authors consider 
a square lattice where, besides the four qubits composing the network, each node has 
one extra qubit (called a processing qubit); see figure [T6J The nodes are divided in three 
categories, black, red, and blue, in such a way that each red or blue node is surrounded 
by four black nodes. After initializing the processing qubits of black nodes in the state 
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Figure 16: Scheme to distribute entanglement in a 2D network based on the surface- 
error correction code. The nodes of the lattice are divided in black, red, and blue nodes. 
The entanglement shared between the nodes is used to measure stabilizer operators. 
Each red (blue) node performs X (Z) measurements on its black neighbouring nodes 
according to the circuit described in figure [I7J 



|0), the entanglement shared between neighbouring nodes is used to perform stabilizer 
measurements. According to the circuit depicted in figure [I7J red nodes perform the 
stabilizer measurements <S>j^ r Xj\f r , while blue nodes measure (£)j^ b Z^ b {N r j, denotes the 
neighbours of red or blue nodes). If the entanglement shared in the network is perfect, 
that is, if the links are given by Bell states, the state of the black processing qubits 
after the stabilizer measurements is transformed into an eigenstate of the surface code. 
The last step of the protocol consists in measuring all black processing qubits except 
the ones that are held by Alice and Bob. The basis for these measurements is chosen to 
create a maximally entangled state between Alice and Bob. 

Due to imperfections, however, errors in the entanglement shared in the network and 
in the operations result in incorrect stabilizer-measurement outcomes. In order to detect 
these errors, the stabilizer measurements are repeated ./V times. For each measurement 
run, a new entangled pair of qubits must be generated between neighbouring sites. 
Since the desired state is an eigenstate of the stabilizer measurements, consecutive 
measurements with differing outcomes indicate an error. Finally, the errors are corrected 
by pairing the error syndromes through the network-based error correction. 

The described protocol tolerates an error rate of approximately 1.67% in the 
quantum channels composing the network, which, in turn, corresponds to links composed 
by Werner states with x ~ 0.98. One of the advantages of the scheme presented in [78] 
is that the same black processing qubit can be used to measure its four neighbouring 
qubits. 

4-2.3. Entanglement distribution with constant resources In the error-correction based 
protocols described above, the number of connections between neighbouring nodes has 
to increase with the size of the network in order to generate long-distance entanglement. 
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Figure 17: Circuit to perform stabilizer measurements using entanglement shared 
between neighbouring nodes. Each wave represents an entangled state. The 
measurement outcomes of the stabilizer measurement (e.g. ® N Xn t ) is given by the 
product of the measurement outcomes at the neighbouring nodes (e.g. x± x x 2 x ^3 x £4). 



On the contrary, it is shown in [TTJ [79] that a constant number of connections is sufficient 
to achieve this task. In the following, we discuss the proposal of [77] which uses cluster 
states in a 3D lattice. We then briefly describe the main idea of [79J, which works both 
for 2D and 3D lattices. 

Three-dimensional cluster states The first protocol achieving long-distance entangle- 
ment with a constant number of connections was proposed in [77]. In that article, it 
was shown that if the fidelity of the connections is larger than a critical value F c , then 
entanglement can be generated between two qubits A and B lying on opposite faces of a 
simple cubic lattice of size N 3 , with N — > 00. In this case, however, it must be stressed 
that the local quantum operations are assumed to be perfect, which is not necessary for 
the other strategies that are fault-tolerant. In what follows, we describe in detail how 
this protocol works, but let us first define the cluster state, which lies at the heart of 
this proposal. 

A cluster state \C) is an instance of graph states [90], and it is usually constructed 
by inserting a qubit |+) at each vertex of the graph and by applying a controlled-phas^j] 
between all neighbours. This state obeys the eigenvalue equation K u \C) — \C) for all 
vertices u, where K u is the stabilizer 

K U = X U Yl Z v , (4.3) 

and where N(u) stands for the neighbourhood of u in the graph. In [77], the desired 
controlled-phases are non-local quantum operations, and therefore they have to be 

\ Under a controlled-phase operation, the state |11) is multiplied by —1, while |10), |01), and |00) 
are unchanged. 
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performed indirectly through the use of the entangled connections and by applying 
some generalized measurements at the nodes. Because the connections are noisy, errors 
are introduced into the system, i.e., the cluster state is not perfect anymore. It can be 
shown that this results in local Z errors occurring independently at the nodes with a 
rate p ~ 6 e. 

Three-dimensional cluster states are known to have an intrinsic capability of error 
correction. In fact, long-range entanglement was shown to be possible between two faces 
of an infinite noisy cubic cluster state [91]. The difference between [91] and [77] is that 
local quantum operations are assumed at every node in the latter case, in particular in 
the two faces. The error correction runs as follows. First, all qubits but A and B are 
measured in either the X or the Z basis. The measurement pattern is such that, in the 
ideal case, the qubits A and B are maximally entangled: 

X A X B \ip AB ) = \ x \4>ab), ^ ^ 

Z A Z B \ip AB ) = \ z \ip AB ), 

where the eigenvalues Ax, £ {— 1, +1} depend on the measurement outcomes x and z 
and are calculated from the stabilizer equations (|4.3p . The effect of the local Z errors is 
to change the sign of these eigenvalues, thus ruining the quantum correlations if no error 
recovery is performed. Then, the lattice is virtually divided into two interlocked cubic 
sublattices (one for each correlation Ax and Xz), and a parity syndrome is assigned to 
nearly all vertices u: 

= n Kv = n Xv n ( 4 - 5 ) 

veAf(u) veM(u) weM'(u) 

where A/ 7 designates the neighbourhood in the corresponding sublattice. Since this 
equation arises from a product of stabilizers, we have that s = 1 if no noise is present 
in the system. However, Z errors do not commute with X measurements, and the 
construction is such that an error changes the sign of two syndromes that are neighbours 
in the other sublattice. 

We are thus back to the error recovery described in section 14.1.31 with two 
differences nonetheless. First, of course, the lattice is three-dimensional and not planar. 
Second, some syndromes lying on the opposite faces cannot be assigned a value because 
of the specific measurement pattern, whereas all syndromes are known in the usual 
error correction. This results in imperfect quantum correlations, but the Monte Carlo 
simulations performed in [77] indicate that the distant qubits A and B are entangled as 
long as the error rate p is smaller than the threshold p c ~ 2.3%. 

Teleportation-based protocol In [73] it is shown that the problem of transmitting 
entanglement in a two-dimensional network is equivalent to the problem of fault-tolerant 
quantum computation in a one-dimensional array of qubits restricted to next-neighbour 
gates. Let us discuss how it works. First, consider a square lattice in which we would 
like to create long-distance entanglement, and let us suppose for the moment that all 
links of the network are perfect. Consider one of the diagonal array of qubits of this 
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(«) , t = l W 




Figure 18: The task of distributing entanglement in a network can be translated to a 
quantum-computation on a line with next-neighbour gates, a) In a fixed time t, the 
qubits composing the computer are given by a diagonal in the lattice. Time evolution 
is translated in a teleportation process between diagonals in the lattice, b) Two-qubit 
gates between times t and t + 1 can be applied by teleporting neighbour qubits to the 
same node in an intermediate diagonal. 

lattice as a one-dimensional quantum computer at the initial time t = 0; see figure fT8h. 
By teleporting each qubit of the computer twice, first right and then above, the state 
of the computer is mapped to an upper-right diagonal array of qubits. The computer 
is now at time t = 1. If between t — to t = 1 one has to implement a two-qubit gate 
between two neighbour qubits, one first teleports one of the qubits up and the other 
right; see figure [18b. In this way, they end up at the same location, where now the gate 
can be locally applied. Finally, each qubit is further teleported to the diagonal defining 
t = 1. 

In the case that each link corresponds to a Werner state, the teleportation scheme 
can be seen as a quantum computation where errors occur. In this way, a fault-tolerant 
quantum computation scheme should be used. In fact, such a scheme exists, where two 
qubits per site are used [ST]. Note that in fault-tolerant error-correction schemes, one 
starts by encoding a known state. However, in practice, in the very first stage the qubits 
are exposed to noise and thus decohere. Furthermore, in order to create an entangled 
state between two distant nodes, one would need a decoding scheme that is also fault- 
tolerant. In [7j5] the authors develop a fault-tolerant encoding- decoding scheme into a ID 
concatenated code jHTJ that works for unknown states. Both the encoding and decoding 
protocols can be done in a one-shot way by measuring syndrome operators (products of 
X and Z operators) and using error correction schemes based one syndrome patterns. 
Thus, by combining the encoding-decoding scheme with the fault tolerant computation 
on a line, one can establish long distance entanglement in the 2D square lattice network. 

A simpler scheme of entanglement distribution in a 3D lattice was also discussed 
in [79], where the authors use a 2D topological code. 
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4-3. Conclusion 

In this chapter, we have seen that long-distance entanglement can be generated in 
mixed-state networks. To this end, some information about where the noise (bit-flip 
and phase errors) enters the system must be collected. This is done by computing a 
series of parity checks at the nodes, which creates a syndrome pattern. Then, the errors 
are corrected by applying local unitaries that are determined by pairing, in an optimal 
way, the detected syndromes. Note that all proposed protocols are based on theoretical 
results for lattices of infinite size, so that it would be interesting, if not necessary, to 
investigate their efficiency for realistic networks of small or medium size. 
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5. Networks with a complex structure 

Up to this point, we have seen protocols that generate entanglement over a large 
distance in quantum networks with a regular structure. In fact, as motivated in 
the Introduction, the creation of remote entangled qubits is of primary importance 
in quantum cryptography, and the first real quantum networks will certainly have 
a lattice structure: Nodes will be spread homogeneously over a plane (the Earth's 
surface, or between buildings in a city) and neighbours will be connected by a quantum 
channel. However, depending on the demand for quantum communication, new links 
will be created between previously unconnected nodes, leading to networks with a richer 
topology. It is therefore a natural task to study the distribution of entanglement in 
complex networks. These networks describe a wide variety of systems in nature and 
society [921 [11], modelling, among other things, chemical reactions in a cell, the spreading 
of diseases in populations, or the connections between routers and computers in the 
Internet. A nice review on this topic can be found in [9]. 

The first approach to networks with a disordered structure was proposed by 
Solomonoff and Rapoport in the early 1950's [93] by combining probabilistic methods 
and graph theory. They introduced the random graph, computing quantities such as 
the critical density and the size of the giant cluster. A decade later Erdos and Renyi 
studied the random graph in depth [911 ESI [96], initiating random graphs as a field of 
mathematics. 

A great deal is known about the random graph, but it lacks many of the important 
features of real-life complex networks. However, it is useful as a starting point, not only 
because it is easier to analyse, but because it exhibits at least two of the most important 
features of complex networks: 1) It possesses the small-world property, which means 
that the length of the shortest path between two nodes increases slowly with the system 
size. 2) It exhibits critical phenomena, with a critical point and a single cluster with 
macroscopic density. A quantum version of the random graph was proposed in [97J, 
where it is shown that its properties change completely when they are subject to the 
laws of quantum physics (section I57T]) . 

Many measures have been suggested to quantify the main properties of the complex 
networks found in the real life, but three concepts seem to occupy a prominent place [9] : 
the small- world [HE], clustering [8], and scale- free [HH] behaviours. A graph with the 
small world property typically has shortest paths between nodes that scale like ln(iV), 
where N is the number of nodes. Graphs with the small-world property may or may 
not have properties such as community structure or some regularity. But they have 
some links, perhaps a small fraction, that connect random nodes, or distant nodes 
if there is community structure. A graph with a large clustering coefficient has the 
property that, if A is connected to both B and C, then it is likely that B is also 
connected to C . A scale free graph has the property that the probability that a 
node has k links decays as a power of k. While the Erdos Renyi random graph 
satisfies the small- world property, it is not scale- free and has zero clustering coefficient. 
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(a) p = 0.10 (b) p = 0.25 (c) p = 0.50 

Figure 19: Evolution process of a random graph of size N = 10: starting from isolated 
nodes, we randomly add edges with increasing probability p, to eventually get a complete 
graph for p = 1. 

Many mathematical models have been introduced over the years to include the two 
other characteristics. However, in all of these networks, as in lattices, the method of 
entanglement percolation is applicable and the percolation thresholds are enhanced by 
some quantum strategies [100] : see section 15731 In the following section, we treat yet 
another, but related, critical phenomenon: the appearance of connected structures of a 
given shape in a random graph. 

5.1. Random graphs 

Here, we briefly introduce random graph theory; the interested reader is referred to [H] 
(and references therein) for a more detailed description of these graphs and a rigorous 
discussion of their properties. 

The theory of random graphs considers graphs in which each pair of nodes % and 
j is joined by a link with probability pij. In the simplest and most studied model, the 
probability is independent of the nodes with p^j = p. The generated graph, which we 
call the Erdos-Renyi (ER) model, is denoted Gn, p and can considered to be the result 
of an evolution process: starting from N isolated nodes, random edges are successively 
added and the obtained graphs correspond to a larger and larger connection probability; 
see figure 

5.1.1. Appearance of subgraphs One of the main goals of random-graph theory is to 
determine the probability p at which a specific property of a graph Gn, p typically arises, 
as N tends to infinity. Many properties of interest appear suddenly, i.e., there exists a 
critical probability p c (N) such that almost every graph has this property if p > p c {N) 
and fails to have it otherwise. Such graphs are said to be typical. For instance, it 
was shown that G^, P is fragmented into small isolated clusters if p < N" 1 , whereas 
percolation occurs above this threshold, that is, one single giant component forms in 
the network. 

A subgraph F = (V, E) of Gn, p is defined as a collection of n < N vertices connected 
by I edges. Like the giant cluster, the subgraphs have distinct thresholds at which they 
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Table 1: Some critical probabilities, according to (15. ip . at which a subgraph F appears 
in random graphs of N nodes connected with probability p ~ N z . For instance, simple 
subgraphs appear at a small connection probability, whereas cycles and trees of all orders 
emerge at the critical value z — — 1. After Albert and Barabasi [9]. 



typically form. It is proven in |101] that the critical probability p c for the emergence of 
F is 

p c (N) = cN~ n l\ (5.1) 

where c is independent of N. It is instructive to look at the appearance of subgraphs 
assuming that p(N) scales as N z , with z G (—oo, 0] a tunable parameter: as z increases, 
more and more complex subgraphs appear; see table [TJ In particular, only node-to-node 
connections appear in the regime z = —2, whereas complete subgraphs (of order four or 
more) emerge above the percolation threshold z — — 1. 

5.2. Quantum random graphs 

A natural extension of the previous scenario to a quantum context was considered in [97J . 
For each pair of nodes, the probability pij is replaced by a quantum state pij of two 
qubits, one at each node. Hence, every node possesses N — 1 qubits that are pairwise 
entangled with the qubits of the other nodes. The study is restricted to the pure-state 
scenario and the pairs of particles are identically connected, so that pij = \<p). A 
quantum random graph is then defined as 

N 

\Gn, p ) = (g) \<p)iy (5.2) 
i<j=l 

Expanding all terms of this expression in the computational basis, one notes that this 
state is the coherent superposition of all possible simple graphs on N nodes, weighted 
by the number of states |11) they possess. For instance, drawing a line for the state 
1 11)0 and nothing for |00), the quantum random graph on three nodes reads: 

\g 3 , p ) = ^ \ @ \) + ^ ( | o \) + \ZJ + l/o)) 

+v^i ( I A) + I A) + I A)) + v^T 3 1 A)- (s-s) 

The analogy with the classical case is that one lets the degree of entanglement of 
the connections scale with the number of nodes: E(ip) = 2ipi ~ N z . In this case, the 
"classical" strategy is to optimally convert each connection of the graph, individually, 

f Note that, contrary to the rest of the Review, a line denotes here a separable state rather than 
entangled qubits. 
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into the Bell pair |$ + ), as described in section l2.2.2l The task of determining the type of 
maximally entangled states remaining after these conversions is mapped to the classical 
problem, and one obtains again the results of table Q] with p = E(tp). 

5.2.1. A complete collapse of the critical exponents The main result in [97J states that, 
in an infinitely large quantum random graph and in the regime p ~ N~ 2 , the state 

i 



which corresponds to any subgraph F = (V, E) composed of n vertices and / edges, can 
be generated with a strictly positive probability using LOCC only. All critical exponents 
of table[T]thus collapse onto the smallest non-trivial value z = —2. This clearly indicates 
that the properties of disordered graphs change completely when they are governed by 
the laws of quantum physics. It is not the purpose of this Review to prove this result, 
but let us describe the main quantum operation that lies behind its proof. In fact, it 
is another example of a generalized measurement on many qubits that enhances the 
distribution of entanglement in quantum networks. 

A joint measurement at the nodes As for multipartite entanglement percolation 
(section 13.2.31) . allowing strategies that entangle the qubits within the nodes yields 
better results than a simple conversion of the links into Bell pairs. In the current case, 
the construction of a quantum subgraph \F) is based on an incomplete measurement of 
all qubits at each node. More precisely, the measurement operators are projectors onto 
the subspaces that consist of exactly m states |1) out of M = N — 1 qubits: 



where n m denotes a permutation of the qubits. Applied on all nodes of the quantum 
random graph, this measurement generates a highly entangled state shared by a random 
subset of nodes. The quantum correlations corresponding to \F) are then extracted by 
a suitable series of local operations at these nodes [97]. Note that this measurement 
allows not only the creation of any quantum subgraph, but also the generation of some 
important multipartite states, such as the GHZ states |43j . 

5.3. Percolation 

One may ask what pure-state entanglement percolation looks like on complex networks 
rather than regular lattices. It is obvious that CEP on complex networks works exactly 
as it does on regular lattices; but it is not obvious how to design pre-processing for 
QEP. However, a particular QEP protocol has been applied with success to a variety of 
complex networks with double-bond, pure, partially-entangled states |lU0j 1102] . As in 
the many cases of links with multiple pairs that we have seen, this setup makes possible 



\F) = ®\$ + ) Ei , 



(5.4) 




(5.5) 
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Figure 20: g-swap protocol. The red dots and lines are involved in one of five swaps that 
transform star into the cycle. Solid lines are partially entangled pure states. Dotted 
lines are post-swapping states. 

the systematic application of a simple local transformation, with the hope that general 
global effects can be understood. 

Following [102j . we consider here networks where each link between nodes consists 
of two pairs in the state \(p); see (12. 5p . CEP consists of the optimal conversion of each 
pair of pairs to a singlet, followed by swapping. The QEP protocol consists of optionally 
performing g/2 swaps at a node of degree g, with q even. The swaps join pairs of links 
terminating at the node. Thus, the node is replaced by a cycle whose bonds are occupied 
with states according to the random outcome of the swapping; see figure [201 We call 
the original node and its bonds a g-star, the transformed object a g-cycle, and the 
transformation a g-swap. Whether it is advantageous to apply the g-swap depends on 
the details of the network. 

The interesting parameter regime for complex networks is often the tree- like regime- 
that is, when there are so few links that the probability that there is more than one path 
connecting two nodes is negligible. For instance, the percolation transition of the ER 
network is in the tree-like regime. The usual way to analyse tree-like networks is with 
generating functions and recursion, which gives exact results in this limit. Although 
the g-swap creates loops, they are small and isolated, so that an entire g-cycle can be 
treated at one step in the recursion. Still, one has a choice in applying g-swap because, 
for instance, it cannot be applied at neighbouring nodes. In |102] this calculation was 
done for a breadth-first application, that is, applying the g-swap at neighbour of a 
starting node, and then at next-nearest neighbour's, where possible, etc. Exact results 
were obtained for various networks. For instance, it was shown that any application 
of g-swap on a tree lowers the threshold with respect to CEP; see section 13.1.11 The 
generating-function analysis shows that, in the tree-like regime, the optimal strategy 
with node-degree k is to either always or never perform g-swap depending on the lattice 
and k. For other lattices this fact comes into play. For example, on the ER network, 
the optimal application of g-swap yields the best results for average degree near 4 with 
a 20% reduction in the threshold, with similar results for the un-correlated scale-free 
network. 

The presence of small cycles is recognized as a key feature of many complex 
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networks. In a network of scientific collaborators, for instance, one can identify small 
groups, each member of which has been a coauthor with each of the others. Of course, 
this phenomenon cannot be modelled in the strict tree-like limit. Thus other models 
and techniques of analysis, largely numeric, are required. For instance, the Watts- 
Strogatz model begins with a one-dimensional chain of nodes connected with nearest- 
neighbour's and next-nearest-neighbour connections. The terminal node of each bond 
is then randomly rewired with small probability p. The original local connections give 
the model community structure, while the rewired bonds provide a few long-range links 
that drastically reduce the shortest path between nodes. Numerical simulations showed 
less impressive improvements than the tree-like models. Still, the critical threshold is 
lowered by g-swapping. However, this form of QEP is not always advantageous. For 
large enough initial entanglement per link, the CEP produces a singlet with probability 
1 and the giant connected component is larger for CEP than for g-swap. 

5.4- Mixed state distribution 

Here we discuss distributing full-rank mixed states on complex networks. One approach 
is to admit that the exponential decay of entanglement due to swapping imposes an 
upper limit on the number swaps that may be performed, before all entanglement 
is lost. The distance corresponding to this limit can be included in calculations of 
statistical properties of entanglement. For instance, if this limit is smaller than the 
correlation length, then the network is essentially fractured with respect to entanglement 
via direct swapping. This approach was taken numerically and with generating functions 
in [102] . On the other hand, it is a priori possible to approximate maximally entangled 
links between distant nodes by concentrating links from ever larger numbers of paths. 
Studies to date have examined possible building blocks to this end |1U3] . This work 
quantifies the gain in concurrence obtained between nodes by employing concentration 
and distribution protocols on multiple paths. The most detailed results were obtained 
for the simplest case of the single purification protocol (SPP), in which the shortest 
path Vab between A and B is identified, and then a shortest path between two nodes 
on Vab is found. Both paths can be used to achieve a final entanglement between A 
and B that is larger than that achieved by using Vab only; see figure [2TJ 

In particular, we examine here the results of the SPP protocol applied to the ER 
random graph |103j . The initial network is the same as in the previous section, except 
that the links are Werner states (I2.8P rather than pure states. As mentioned above, it 
is impossible to extract a pure state from a finite number of Werner states. We instead 
search for a protocol that produces the mostly highly entangled mixed states possible. 
More specifically, we seek to maximize the average of the concurrence ( 12. lip . 



where the sum is over all pairs of nodes in the network, it is the probability that the 
protocol connecting a and (3 was successful and C{x a ^) is the resulting concurrence 




(5.6) 
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Figure 21: Establishing entanglement between nodes A and B. (a) The shortest path 
Tab between A and B; the geometry of the path is irrelevant, so we represent it by a 
straight line with individual links not shown. Other paths connecting A and B are not 
shown, (b) The shortest path Vab (solid line with a dashed segment) between A and 
B. Between the endpoints of subpath S (dashed segment) there is an alternate path A 
(dotted line). 
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Figure 22: Regions in which the single purification protocol is advantageous v.s. rescaled 
Werner parameter y and path length a, and for various values of excess alternate path 
length b. The protocol yields higher average concurrence inside the closed curves. From 
the outermost to innermost curve the values of b are 0, 0.01, 0.07, 0.11, 0.135. The curve 
cutting through the closed curves is y a = 2y and maximizes (independently of b) the 
increase in average concurrence with respect to a. 

between the pair. 

The most naive method to entangle two nodes A and B is to perform entanglement 
swapping repeatedly between Werner states along the shortest path Vab joining A and 
B according to (12. 18ft : see figure l2"Th. However, in some cases, higher entanglement may 
be obtained by additionally concentrating the entanglement from an alternate path 
connecting intermediate nodes on the path Vab as shown in figure [2Tb . In this single- 
purification protocol, one swaps along the subpath S as well as the shortest available 
alternate path A, then performs a purification on the two resulting parallel states, and 
finally performs swapping at all remaining links. Whether a higher entanglement results 
on average than that from simply swapping along Vab depends on the lengths of the 
paths and the Werner parameter x. 

For instance, consider the parameters L = \ \Vab\\, a = \\S\\/L, b = (||^4|| — ||«S||)/L, 
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Figure 23: Regions in which a serial multiple purification protocol is advantageous v.s. 
rescaled Werner parameter y and total fractional path length a, and for various numbers 
of subpaths n. The multiple path protocol is advantageous only between the solid and 
dotted lines. 



and y = x x l L . A detailed analysis of the region in this parameter space for which this 
purification protocol yields a higher average concurrence is given in [103] , and is shown 
graphically in figure [22J 

Of course a network may offer more possibilities than concentrating the 
entanglement from a single neighbouring path. One could repeatedly concentrate 
entanglement from paths in parallel or series. Let us consider the latter case, where 
instead of a pair of paths S and A, we have n pairs and that together the n paths Si 
cover a fraction a of Vab- For simplicity we take ||<S|| = ||-4.||. The region in parameter 
space for which this protocol is advantageous were computed exactly and are shown 
in figure [221 One sees that many short paths are better than few long paths, but that 
the maximum extent of the good region approaches a limit as n — > oo. 

Suppose we apply only the single-purification protocol at every possibility on a 
network. Computing the average concurrence over a random network is in general 
difficult. However, the ER network at its critical point Np = 1 is in the tree-like regime 
and has further statistical properties that facilitate calculations. In particular, the exact 
asymptotic increase in concurrence over the naive swapping protocol as a function of 
Werner parameter x and number of nodes iV was found to be AC ~ AN~ 2 (1 — x)~ 4 . 
Here A w 6.5 x 10 -5 is a constant that is easily computed numerically. Although this 
expression diverges as x — >■ 1, this occurs outside the asymptotic regime where the 
purified paths have lengths much shorter than the radius of the giant cluster. 

It is important to note here that the above protocol is not capable of entangling 
arbitrarily distant nodes independently of the initial concurrence per link, but it is one 
that can be applied in principle to arbitrary networks. 

The preceding protocols, although treating mixed states, have assumed perfect 
operations. If we further assume the noise model described in (12.211) and (12.241) we find 
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that the protocol is rather sensitive to noise. The results are quantitatively similar to [6] 
with the advantage destroyed for noise levels larger than a few percent. For instance, 
the maximum average concurrence gain from the single purification protocol AC was 
found to be 



The richer variety and topologies of complex networks provide a fertile terrain for 
developing more sophisticated distribution protocols. The studies reviewed here make 
only the first steps. In particular, they apply relatively simple quantum procedures for 
which averages over quantum outcomes and classical disorder may be performed. Other 
important directions remain completely unexplored, such as the dynamic creation and 
distribution of entanglement. 




(5.7) 
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6. Conclusion 

In this review we have given an account of theoretical progress on the distribution of 
entanglement in quantum networks. We have seen that this inquiry has been driven 
by the experimental results on the building blocks of entanglement distribution. These 
results are impressive and promising, but also make evident fundamental and technical 
barriers. We have also seen that the task of distribution is intimately connected with 
theoretical questions about the nature of entanglement. In this setting, these questions 
are focused on the extent to which it can be measured and inter-converted. While the 
basic concepts (direct transmission, swapping, purification, error correction) have been 
present for over twenty years, their application to real networks is still in an exploratory 
phase. 

The firsts steps to transmission of entanglement embodied in a variety of physical 
systems have already been taken. In fact, the most important application to date, the 
distribution of quantum cryptographic keys, has been demonstrated. Teleportation, as 
well, has been demonstrated in several systems. But the problem of exponential decay 
of fidelity has not been solved. Methods using the geometry of higher-dimensional 
networks to effectively concentrate entanglement from some sections for use in others 
shows promise in overcoming this difficulty At the same time, they show strong 
connections with the theory of classical networks and graphs. Classical error correction 
also has proven to be useful in this regard, with interesting quantum connections, such 
as defining syndromes via weak measurements that preserve information. 

It is clear that the study of entanglement distribution on mixed-state networks will 
be of prime importance. For instance, quantum cryptography will continue to be one 
of the main technologies driving research. Quantum repeaters are expected to become 
more robust, thereby allowing small networks to form. This will require new protocols 
adapted to these small networks. Eventually, we expect to see application of pioneering 
work in entanglement percolation and error correction. Finally, as the number of nodes 
increases to the point that statistical methods can be applied, we expect to see vigorous 
activity in the theory of complex quantum networks. 
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